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The  objective  of  the  present  program  was  to  study  the  aerothermochemical 
structure  of  laminar  premixed  and  nonpremixed  flames  through  (a)  non-intrusive 
experimental  determination  in  reduced  and  elevated  pressure  environments,  (b) 
computational  simulation  using  detailed  flame  and  kinetic  codes,  and  (c) 
asymptotic  analysis  with  simplified  and  reduced  mechanisms.  Useful  theoretical 
and  experimental  contributions  were  made  on  the  determination  of  the  burning 
rates  and  flame  kinetics  of  the  lower  hydrocarbons,  on  the  understanding  of  the 
physical  and  chemical  parameters  influencing  soot  formation  in  diffusion 
flames,  on  the  identification  of  the  role  of  kinetics  and  system  non- 
adiabaticity  in  flammability  limits,  and  on  adiabatic  flame  stabilization. 

These  results  are  relevant  to  the  fundamental  and  practical  issues  of  flame 
kinetics,  turbulent  combustion,  soot  formation,  radiative  heat  transfer, 
flame  extinction,  stabilization  and  flammability,  and  supersonic  combustion. 
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SUMMARY/OVERVIEW 

The  objective  of  the  present  program  is  to  study  the 
aerothermochemical  structure  of  laminar  premixed  and  nonpremixed 
flames  through  (a)  non-intrusive  experimental  determination  in  reduced  and 
elevated  pressure  environments,  (b)  computational  simulation  using  detailed 
flame  and  kinetic  codes,  and  (c)  asymptotic  analysis  with  simplified  and 
reduced  mechanisms.  Useful  theoretical  and  experimental  contributions 
have  been  made  on  the  determination  of  the  flame  burning  rates  and  flame 
kinetics  of  the  lower  hydrocarbons,  on  the  understanding  of  the  physical 
and  chemical  parameters  influencing  soot  formation  in  diffusion  flames,  on 
the  identification  of  the  role  of  kinetics  and  system  nonadiabaticity  in 
flammability  limits,  and  on  adiabatic  flame  stabilization.  These  results  are 
expected  to  be  useful  in  the  general  interest  of  AFOSR  in  the  fundamental 
and  practical  issues  of  flame  kinetics,  turbulent  combustion,  soot  formation, 
radiative  heat  transfer,  flame  extinction,  stabilization  and  flammability,  and 
supersonic  combustion. 


ACCMPLISHMENTS 

The  highlights  of  our  accomplishments  can  be  found  in  the  annual 
reports  submitted  to  the  program  director,  as  well  as  the  journal  papers 
which  have  either  appeared  or  have  been  accepted  for  publication.  Thus  only 
a  very  brief  summary  of  those  completed  works  of  archival  value  are 
discussed  in  the  following. 

1.  Determination  of  Burning  Velocities  and  Chemical  Kinetics  of 
Hydrocarbon/Air  Mixtures 

A  major  accomplishment  of  our  AFOSR  program  has  been  the 
development  of  the  counterflow  twin  flame  technique  which  allows  the 
accurate  determination  of  the  laminar  flame  speeds  of  combustible  mixtures. 
Basically,  practically  all  existing  techniques  for  the  determination  of  laminar 
flame  speeds  impart  stretch  and/or  heat  loss  effects  on  the  data  obtained, 
causing  significant  and  systematic  spreads  in  their  values.  Thus  when  such 
flame  speed  data  are  used  to  calibrate  or  validate  chemical  kinetic 
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mechanisms  and  rate  constants,  the  fundamental  chemical  information  so 
determined  can  be  significantly  “contaminated"  or  even  falsified.  In  the 
technique  we  proposed,  such  effects  can  be  systematically  subtracted  out, 
yielding  flame  speed  data  of  unusual  fidelity. 

For  the  present  grant  period,  we  have  capitalized  on  our  expertise  in 
this  area  in  the  study  of  a  variety  of  fuel/oxidizer  systems.  First,  we  have 
extended  our  previous  study  of  methane/air  flames  to  cover  all  the  C2 
hydrocarbon  fuels  of  ethane,  ethylene,  and  acetylene,  as  well  as  propane. 
Such  an  extension  is  crucial  because  methane  only  has  the  C-H  bonds 
although  the  breaking  of  the  C2  bonds  is  important  in  the  kinetics  of  all  of 
the  higher  hydrocarbons.  Through  this  study  we  have  found  that  all  existing 
C2  kinetic  schemes  are  inadequate  to  describe  our  experimental  results. 
Consequently  a  new  kinetic  scheme  which  was  able  to  correlate  practically 
all  of  our  experimental  results  was  proposed.  It  is  also  important  to  mention 
that,  through  our  effort  over  the  years,  there  now  exist  reliable  experimental 
flame  speed  data  for  methane,  ethane,  ethylene,  acetylene,  and  propane, 
determined  over  extensive  ranges  of  stoichiometry  and  up  to  pressures  of 
about  3  atmospheres.  These  data  are  essential  for  both  fundamental  and 
practical  studies.  This  work  is  reported  in  Publication  No.l. 

We  have  also  studied  the  oxidation  kinetics  of  hydrogen.  The  practical 
interest  to  AFOSR  is  on  supersonic  combustion,  while  the  fundamental 
interest  is  that  hydrogen  oxidation  forms  a  fundamental  building  block  for 
the  oxidation  of  carbon  dioxide  as  well  as  the  hydrocarbons.  In  this  study  we 
have  focused  on  very  lean  mixtures  because  flame  speed  data  for  these 
situations  do  not  exist,  and  because  of  the  interest  in  low-  and  intermediate- 
temperature  kinetics  which  control  these  weakly-burning  flames.  In  may  be 
noted  that  previous  studies  on  very  lean  hydrogen/air  flames  have  all  been 
complicated  by  the  presence  of  flamefront  instabilities  which  are  favored  to 
occur  for  these  high  Lewis  number  flames.  However,  since  such  instabilities 


are  suppressed  by  the  positively  strained  flow  field  of  the  counterflow  flame,  p 
we  have  succeeded  in  obtaining  data  for  smooth  hydrogen/air  flames.  Our  jjp 
experimental  results  show  that  there  is  something  seriously  wrong  with  the 
existing  hydrogen /oxygen  kinetics,  in  these  low  flame  temperature  regimes,  t — - — » 
in  that  the  predicted  flame  response  obtained  by  using  existing  mechanisms 
completely  fails  to  predict  the  observed  behavior  even  qualitatively.  It  seems  } 
that  the  actual  kinetics  is  more  robust  than  indicated.  This  could  be  caused  * 
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by  the  great  uncertainty  associated  with  the  kinetics  involving  the  HO2  and 
H2O2  radicals,  as  generally  acknowledged  by  the  chemical  kinetics 
community.  This  work  is  reported  in  Publication  No.2. 

In  an  effort  to  further  scrutinized  the  kinetics  of  HO2  and  H2O2.  we 
have  performed  additional  experiments  on  methanol  oxidation,  for  which 
these  radicals  are  expected  to  exert  different  influences.  Some  useful 
insights  have  been  gained,  although  the  issue  of  hydrogen  oxidation  remains. 
The  work  on  methanol  oxidation  is  reported  in  Publication  No. 3 

The  laminar  flame  speed  data  we  have  been  reporting  over  the  past 
few  years  have  now  enjoyed  wide  popularity  among  combustion  kineticists 
and  flame  modellers.  Since  these  data  are  scattered  in  various  journal 
papers  and  also  reported  in  different  formats,  upon  request  from  many  of 
our  colleagues  we  have  consolidated  and  re-plotted  all  our  data  in  a  uniform 
format.  The  compilation  is  reported  in  Publication  No.  4,  which  will  appear 
in  a  monograph  on  reduced  mechanisms.  All  contributors  to  this 
monograph  are  asked  to  compare  their  calculated  results  against  our 
reported  flame  speeds. 

We  next  note  that  in  efforts  to  describe  the  effects  of  chemical 
reaction  on  combustion  phenomena  with  greater  realism  than  the  one-step 
reaction  frequently  used  in  theoretical  analysis,  various  reduced  mechanisms 
based  on  steady-state  and  partial  equilibrium  concepts  have  been  developed. 
Most  of  these  mechanisms  consist  of  four  steps.  We  have,  however,  found 
that  these  four  step  mechanisms  are  not  adequate  to  describe  the  computed 
flame  responses  using  detailed  chemical  mechanisms.  The  deficiency  was 
identified  to  be  the  steady-state  assumption  for  the  CH3  radical.  A  five-step 
mechanism  was  subsequently  proposed,  and  the  agreement  was  significantly 
improved.  This  result  is  expected  to  be  of  particular  relevance  to 
phenomena  that  are  strongly  dependent  on  accurate  descriptions  of  the  CH3 
such  as  the  formation  of  prompt  NOx  and  fuel  rich  combustion.  This  work  is 
reported  in  Publication  No.  5. 

In  another  work  of  great  interest,  by  using  the  laminar  flame  speed 
data  for  methane/air  mixtures  we  have  extracted  the  kinetic  parameters  of 
the  reaction  order  (n)  and  activation  energy  (E)  for  the  equivalent  one-step 
overall  reaction.  The  results  show  that  these  values  are  far  from  being 
constants.  Instead  they  vary  significantly  not  only  with  the  equivalence  ratio, 
but  also  with  the  system  pressure.  For  example,  the  activation  energy  is 


found  to  continuously  increase  with  pressure.  The  reaction  order  is  found  to 
be  just  about  one  at  atmospheric  pressure  and  decreases  with  increasing 
pressures.  For  highly-diluted  environments  n  can  even  assume  negative 
values.  Such  behavior  can  be  explained  on  the  basis  of  the  competition 
between  termination  and  branching  reactions.  This  work  is  reported  in 
Publication  No.  6. 

The  above  work  are  all  concerned  with  premixed  flames.  We 
have  in  addition  also  performed  extinction  experiments  on  methane/air 
diffusion  flames.  The  extinction  stretch  rates  were  consequently  determined 
as  a  function  of  local  stretch  measured  by  LDV.  The  data  were  then 
compared  with  numerical  simulations  and  useful  insight  into  the  extinction 
mechanism  were  gained.  This  work  is  reported  in  Publication  No.  5. 

We  close  this  section  by  emphasizing  the  significance  of  our 
contributions  in  the  area  of  fundamental  flame  and  kinetics  research.  At  the 
23rd  International  Combustion  Symposium,  the  principal  investigator  was 
awarded  the  Combustion  Institute  Silver  Medal  for  the  development  of  the 
counterflow  twin-flame  technique  for  the  determination  of  the  laminar  flame 
speed. 

2.  Theory  of  Flammability  Limits 

Although  the  term  “flammability  limit"  has  been  extensively  used  in 
the  combustion  literature,  its  real  meaning  has  remained  unclear  and 
undefined.  Thus  a  major  and  sustained  undertaking  of  the  present  program 
is  to  understand  the  phenomenon  of  flammability  limits  and  to  be  able  to 
predict  their  values  from  first  principles.  Through  various  considerations, 
we  have  argued  that  a  truly  fundamental  definition  of  flammability  limits 
should  be  related  to  the  state  at  which  the  classical  one-dimensional  planar 
flame  fails  to  propagate,  as  originally  adopted  by  Spalding  in  his  attempt  to 
define  such  a  limit.  By  introducing  a  volumetric  heat  loss  rate  into  the 
energy  equation,  and  by  assuming  a  one-step  overall  reaction  as  well  as 
conventional  transport  properties,  subsequent  analytic  solutions  yield  an 
extinction  state  in  terms  of  a  turning  point  response  of  the  mass  burning 
rate  to  increasing  heat  loss,  indicating  the  incipient  loss  of  balance  between 
heat  loss  and  chemical  heat  generation.  It  is  further  shown  that  the  ratio  of 
the  limit  flame  speed  to  the  adiabatic  flame  speed  is  e1/2=0.61. 
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The  loss  theory  is  inherently  incapable  of  explaining  the  chemical 
kinetic  dependence  of  the  flammability  limits,  in  terms  of  the  frequently- 
present  chain  branching  and  termination  mechanisms,  on  such  system 
parameters  as  pressure,  stoichiometry,  and  fuel  type.  In  Publication  No.  8 
we  first  proposed  a  chain-based  approach,  in  which  it  was  argued  that 
competition  between  chain  branching  and  termination  reactions  must  be 
important  to  the  weakly- reactive  flames  close  to  the  flammability  limits. 
Thus  by  numerically  simulating  the  one-dimensional  flame  propagation  with 
detailed  chemistry  and  transport  properties,  but  without  radiative  heat  loss, 
a  normalized  sensitivity  coefficient  of  the  termination  reaction  to  variation*3 
in  the  branching  reaction  can  be  calculated.  This  coefficient,  termed  the 
flammability  exponent,  was  found  to  assume  a  value  close  to  unity  in  the 
vicinity  of  the  experimentally-observed  flammability  limits  of  a  great  variety 
of  mixtures. 

It  is  clear  that  the  heat  loss  theory  of  Spalding  and  our  chain 
branching- termination  criterion  are  based  on  completely  different  concepts 
and  approaches,  and  that  neither  one  is  complete.  These  two  concepts  have 
since  been  unified  in  Publication  No.  9.  The  approach  towards  such  a 
unified  treatment,  with  quantitative  accuracy,  is  to  numerically  simulate  the 
freely-propagating,  one-dimensional  flame  with  detailed  chemistry  and 
diffusive  properties,  allowing  for  radiative  heat  loss  from  the  major  species, 
and  to  generate  the  characteristic  turning  point.  At  the  same  time,  the 
relative  efficiencies  of  the  key  termination  versus  branching  reactions  will 
be  assessed  through  the  concept  of  the  flammability  exponent,  especially 
near  and  at  the  state  of  the  turning  point. 

Calculations  have  been  conducted  for  lean  methane/air  and  rich 
hydrogen/air  flames  with  and  without  radiative  heat  loss.  The  results  show 
that  the  nonadiabatic  flame  indeed  exhibits  a  turning  point.  At  the  turning 
point,  the  flame  speed  and  temperature  are  very  close  to  those  predicted 
from  asymptotic  theories  when  referenced  to  the  adiabatic  values,  which  are 
well  established.  It  is  further  shown  that  the  flammability  exponent  assumes 
a  value  close  to  unity  at  the  turning  point. 

The  above  results  provide  a  unified  interpretation  of  flammability 
limits  by  recognizing  that  the  influence  of  heat  loss  and  chain  termination 
on  flame  propagation  become  critical  almost  simultaneously.  As  such, 
flammability  limits  determined  through  the  turning  point  behavior  are 
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almost  identical  to  those  determined  through  the  unity  flammability 
exponent  criterion.  The  result  is  entirely  physically  reasonable  by 
considering  the  following.  As  the  flammability  limit  is  approached,  the 
adiabatic  flame  temperature  is  reduced  due  to  the  change  in  stoichiometry. 
This  weakens  the  temperature  sensitive  two-body  branching  reaction  but 
usually  has  a  much  weaker  influence  on  the  termination  reactions,  especially 
when  it  involves  three  bodies.  Therefore,  a  state  will  be  reached  at  which 
the  termination  reaction  becomes  overwhelmingly  strong,  as  indicated  by 
the  unity  flammability  exponent,  and  consequently  causes  a  rapid  slow  down 
of  the  overall  reaction  rate  and  thereby  the  overall  heat  release  rate.  This,  in 
turn,  makes  the  influence  of  heat  loss  relatively  severe,  and  therefore 
simultaneously  results  in  the  extinction  turning  point. 

3.  Soot  Formation  in  Diffusion  Flames 

We  have  completed  two  projects  related  to  soot  formation  in  diffusion 
flames.  In  the  first  project,  we  have  continued  our  study  of  the  influence  of 
dilution  on  soot  formation.  As  background,  we  note  that  a  major  concept 
developed  in  recent  studies  on  soot  formation  in  diffusion  flames  is  that  the 
flame  temperature  exerts  a  dominant  influence  on  soot  formation.  This 
concept  has  been  experimentally  investigated  by  lowering  the  flame 
temperature  by  diluting  the  fuel  stream  with  an  inert.  We  have,  however, 
argued  that  since  the  act  of  dilution  not  only  lowers  the  flame  temperature 
but  also  its  concentration,  this  concentration  modification  should  also  have 
an  effect  on  soot  formation.  By  subsequently  conducting  the  experiment  by 
holding  either  the  flame  temperature  or  the  fuel  concentration  fixed,  we 
have  shown  that  not  only  temperature  and  concentration  are  equally 
important  in  soot  formation,  in  many  instances  dilution  can  be  even  more 
important. 

Our  previous  experiments  were  conducted  by  using  the  counterflow 
flame,  although  most  of  the  literature  experiments  were  performed  with  the 
co-flow  flame.  We  have  therefore  repeated  our  experimental  investigation 
for  the  co-flow  flame,  using  the  same  isolation  technique  as  developed 
previously.  The  results  confirm  our  previous  finding  of  the  simultaneous 
importance  of  temperature  and  concentration  in  soot  formation.  This  work 
is  reported  in  Publication  No.  10. 
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Our  second  project  on  soot  formation  involves  the  influence  of 
additives.  Carbon  monoxide  and  oxygen  were  added  to  either  the  fuel  or  the 
oxidizer  side  of  a  diffusion  flame,  and  the  chemical,  dilution,  and 
temperature  effects  were  isolated  and  studied.  Results  show  that  the 
addition  of  carbon  dioxide  is  usually  suppressive,  while  the  addition  of 
oxygen  can  be  either  promoting  or  suppressive.  This  work  is  reported  in 
Publication  No.  11. 

4.  Adiabatic  Flame  Stabilization 

The  mechanism  with  which  a  Bunsen  flame  is  stabilized  at  the  burner 
rim  is  considered  to  be  well  established.  The  concept  of  stabilization  is 
based  on  the  existence  of  a  dynamic  balance  between  the  local  flow  velocity 
and  flame  velocity  at  a  certain  point  on  the  flame  surface,  and  the  ability  of 
the  flame  to  adjust  its  flame  velocity,  and  thereby  the  location  of 
stabilization,  through  heat  loss  to  the  rim.  Blowoff  occurs  when  such  a 
balance  can  not  be  achieved  everywhere  on  the  flame  surface.  This 
mechanism  has  also  served  as  the  fundamental  concept  in  flame  stabilization 
through  heat  loss  to  the  stabilization  body  in  other  practical  situations. 

However,  recent  studies  on  the  dynamics  of  stretched  flames  have 
shown  that,  in  addition  to  heat  loss,  the  flame  speed  can  also  be  modified  by 
flow  nonuniformity,  flame  curvature,  and  mixture  preferential  diffusion.  It 
is.  therefore,  of  interest  to  explore  the  possibility  that  flame  stabilization  can 
be  achieved  in  the  absence  of  heat  loss  due  to  these  additional  mechanisms. 
The  configuration  adopted  for  study  is  the  inverted  flame  because  of  its 
symmetry.  By  measuring  the  temperature  distribution  in  the  stabilization 
region,  we  conclusively  demonstrated  that  the  flame  can  indeed  be 
stabilized  in  the  absence  of  heat  loss.  We  subsequently  formulated  a  theory 
that  satisfactorily  describes  the  stabilization  of  this  adiabatic  flame  through 
the  combined  influence  of  flow  nonuniformity  and  flame  curvature. 

It  may  be  emphasized  that  the  concept  and  methodology  developed  in 
this  study  will  be  useful  in  future  studies  of  flame  stabilization  in  combustors, 
as  well  as  the  quenching  of  laminar  flamelets  constituting  turbulent  flames. 

This  work  is  reported  in  Publication  No.  12. 
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5.  Review  Activities 

A  review  article  on  subsonic  and  supersonic  flame  stabilization,  listed 
as  Publication  No.  13,  was  prepared  and  presented  at  a  workshop  on 
supersonic  combustion  at  NASA  Langley  Research  Center. 
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Abstract — An  experimental  and  numerical  study  of  methanol  kinetics  has  been  conducted.  A  detailed 
kinetic  scheme  was  compiled  which  closely  predicts  properties  of  mixtures  of  methanol,  oxygen,  and  inert 
for  a  variety  of  experimental  configurations  and  conditions.  The  scheme  incorporates  the  most  recent 
kinetic  information  and  was  tested  agains'  experimental  data  for  the  propagation  speeds  and  structure  of 
laminar  flames  as  well  as  the  species  concentration  evolutions  in  flow  reactors,  static  reactors,  and  shock 
tubes  The  laminar  flame  speeds  of  atmospheric  methanol/air  mixtures  were  determined  using  the  counter¬ 
flow  flame  technique  over  extensive  lean-to-rich  fuel  concentration  ranges  and  for  initial  mixture  tem¬ 
peratures  rangin',  from  318  to  368  K,  while  the  experimental  data  on  the  laminar  flame  structure  and  from 
reactors  and  shock  tubes  were  obtained  from  the  literature.  The  scheme  compiled  herein  includes  the 
detailed  C,.  C-.  and  methanol  submechanisms  and  yields  close  agreement  with  all  of  the  experimental 
methanol/air  laminar  flame  speeds  as  well  as  previously  determined  laminar  flame  speeds  of  mixtures  of 
CH,  and  the  C- -hydrocarbons  with  air.  The  relative  importance  and  influence  of  the  individual  reactions 
on  the  flame  speed  and  reaction  mechanism  were  assessed  with  the  aid  of  sensitivity  and  species  con¬ 
sumption  path  analyses  The  study  also  demonstrates  that  accurate  prediction  of  laminar  flame  speeds  is 
only  a  necessary  but  not  sufficient  condition  for  the  validation  of  the  mechanism,  and  that  results  from  the 
flame  structure  as  well  as  reactors  and  shock  tubes  are  also  needed  for  further  validation.  Both  flow  and 
static  reactor  oxidation  studies  indicate  that  the  reaction  CH,OH  +  OH  ->  Products  may  have  a  slower 
rate  than  that  reported  in  recent  literature 


1  INTRODUCTION 

In  the  search  of  viable  alternate  fuels  for  reduced  pollutant  emissions  and  lessened 
dependence  on  available  petroleum  supplies,  alcohols  and  their  blends  with  conven¬ 
tional  petroleum  stocks  are  among  the  strongest  candidates.  In  particular,  methanol, 
the  simplest  of  the  alcohols,  has  been  the  subject  of  extensive  laboratory  and  road 
testing  studies  as  well  as  legislative  directives  encouraging  its  use. 

There  are  various  advantages  and  disadvantages  in  the  use  of  methanol  as  a 
practical  fuel  (Bowman.  1975;  Aronowitz  et  al„  1979;  Westbrook  and  Dryer,  1979; 
Gulden  1982;  Dove  and  Warnatz.  1983;  Cribb  etal..  1984;  Grotheer  and  Kelm.  1989; 
Norton  and  Dryer.  1989.  1990;  Norton,  1990).  Specifically,  the  major  advantages  are; 
(a)  convenient  storage  in  vehicles  as  a  liquid  fuel  under  normal  temperature  and 
pressure,  (b)  diversified,  petroleum-independent  sources  of  production,  such  as 
natural  gas.  wood,  coal,  and  renewable  biomass  feedstock,  (c)  reduced  propensity  for 
the  formation  of  pollutants  such  as  soot.  CO.  NO, .  and  unburned  hydrocarbons,  and 
(d)  high  octane  number  rating.  Among  the  various  disadvantages  (Norton,  1990)  are: 
(a)  lower  energy  density,  (b)  propensity  to  form  formaldehyde,  (c)  infinite  miscibility 
with  water  and  thereby  tendency  to  dissolve  atmospheric  moisture,  and  (d)  low  vapor 
pressure  leading  to  cold  start  problems  typically  below  15°C. 

Satisfactory  studies  of  methanol  combustion,  at  both  the  fundamental  and  prac¬ 
tical  levels,  require  a  good  understanding  of  the  chemical  kinetics  of  its  oxidation  and 
the  associated  pyrolysis  steps.  Early  attempts  to  model  methanol  oxidation  were 
fraught  with  significant  uncertainty  until  Tsang  (1987)  compiled  an  extensive  chemical 
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data  base  for  methanol-related  reactions,  with  a  rigorous  evaluation  of  all  available 
data.  Based  on  this  data  base.  Norton  and  Dryer  (1989,  1990)  revised  the  previous 
scheme  of  Westbrook  and  Dryer  (1979)  to  analyze  their  recent  flow  reactor  data  as 
well  as  published  data  obtained  from  static  reactors  and  shock  tubes.  The  funda¬ 
mental  understanding  gained  from  these  studies  is  that,  compared  to  the  oxidation  of 
methane,  the  presence  of  the  C-O  bond  in  the  methanol  molecule  significantly  alters 
its  kinetic  behavior  in  that  completely  different  kinetic  paths  are  being  followed  after 
the  initial  fuel  reactions.  As  such,  the  current  extensive  knowledge  of  the  C,  and  C2 
chemistry  developed  for  methane  oxidation  is  intrinsically  inadequate  to  describe 
methanol  oxidation.  It  is,  however,  also  recognized  that  although  significant  advances 
in  the  understanding  of  methanol  kinetics  have  been  achieved  from  these  worthwhile 
studies,  more  validation  studies  over  wider  temperature,  concentration,  and  pressure 
ranges  are  needed,  especially  for  some  of  the  crucial  kinetic  steps  involving  inter¬ 
mediates  such  as  CH;OH,  CH,0,  CH-O.  and  HCO. 

Such  additional  validations  can  be  achieved  through  simulation  of  flame  propa¬ 
gation  and  structure  because  of  the  extensive  variations  in  the  temperature  and  species 
concentrations  across  the  flame.  Methanol  flames  have  been  the  subject  of  a  number 
of  experimental  and  numerical  investigations  (Wiser  and  Hill,  1955;  Henderson  and 
Hill.  1956;  Gibbs  and  Calcote,  1959;  De  Wilde  and  Van  Tiggelen,  1968;  Akrick  et  al ., 
1978;  Westbrook  and  Dryer.  1979.  1980;  Westbrook,  1980;  Hirano  et  al .,  1981; 
Vandooren  and  Van  Tiggelen,  1981 ;  Giilder.  1982;  Metghalchi  and  Keck,  1982;  Dove 
and  Warnatz,  1983;  Anderson  et  al.,  1984;  Olsson  et  al.,  1986,  1987;  Grotheer  and 
Kelm,  1989;  Pauwels  et  al.,  1989,  1990;  Bradley  et  al.,  1991)  and  certain  rates  have 
been  determined  and/or  adjusted  based  on  comparisons  between  experiments  and 
calculations.  The  fidelity  of  such  simulations,  however,  depends  on  the  accuracy  of  the 
experimental  data.  Specifically,  it  is  generally  recognized  (Law,  1989)  that  while  the 
numerical  simulation  employs  the  one-dimensional  planar  configuration,  the  exper¬ 
imental  flame/flow  configurati.  as  can  deviate  significantly  from  such  an  idealized 
prescription  and  consequently  suffer  substantial,  and  frequently  unquantified,  flame 
stretch  effects.  To  demonstrate  this  point,  the  limited  experimental  da  a  on  the 
laminar  flame  speeds.  S®.  of  methanol/air  flames  under  1  atmosphere  pressure  and 
initial  mixture  temperature.  Tu,  of  298  K  (Wiser  and  Hill,  1955;  Gibbs  and  Calcote, 
1959;  Gulder.  1982;  Metghalchi  and  Keck.  1982)  are  compiled  in  Figure  1  together 
with  the  present  experimental  results;  note  that  S®  has  been  conventionally  ter, r  ed  the 
burning  rate  of  the  mixture.  It  is  seen  that  these  data  not  only  depend  on  the  :  oecific 
technique  used,  whether  flame  tube  (Wiser  and  Hill.  1955),  Bunsen  flame  (Gibi  and 
Calcote.  1959),  and  constant  volume  bomb  (Gulder.  1982;  Metghalchi  and  Keck, 
1982).  but  there  is  also  systematic  difference  between  data  determined  by  the  same 
technique  (Gulder,  1982;  Metghalchi  and  Keck.  1982).  The  spread  of  these  data  is 
sufficientlv  wide  as  to  cause  concern  in  the  assessment  of  the  rate  constants  of  various 
crucial  reactions. 

In  view  of  the  above  consiuerations.  the  first  objective  of  the  present  investigation 
is  to  provide  experimental  data  of  improved  accuracy  on  S®  for  atmospheric 
methanol/air  flames  over  extensive  ranges  of  stoichiometry  and  for  several  initial 
mixture  temperatures.  The  technique  used  is  that  of  the  counterflow  twin  flame,  with 
systematic  elimination  or  minimization  of  stretch  effects  (Wu  and  Law,  1985;  Yu 
et  al .  1986;  Egolfopoulos  et  al..  1989;  Law.  1989;  Zhu  et  al..  1989;  Egolfopoulos  et  al., 
1990.  1991 ).  The  availability  of  these  data  is  essential  not  only  for  the  present  kinetic 
stud:ss  but  also  for  simulation  studies  of  more  complex  fundamental  and  practical 
systems. 

The  second  objective  is  to  compare  our  experimental  S®  data  with  the  numerically- 
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FIGURE  1  Experimental  laminar  flame  speeds.  S?(<j>).  for  methanol/air  mixtures  at  Tu  =  298  K  and 
p  =  I  atm.  from  various  sources. 


calculated  values  by  grafting  the  mechanol  submechanism  of  Norton  and  Dryer 
(1989.  1990)  onto  a  previously  developed  C,  kinetic  scheme  (Egolfopoulos  et  al.,  1990, 
1991)  which  closely  predicts  the  laminar  flame  speeds  of  mixtures  of  CH4  and  the 
C; -hydrocarbons  with  0:  and  N:.  These  comparisons  will  be  used  to  assess  the 
validity  of  existing  rate  coefficients  through  sensitivity  and  species  consumption  path 
analyses.  Such  a  comparison  will  also  demonstrate  that  agreement  in  the  laminar 
flame  speed  is  a  necessary  but  not  sufficient  condition  for  the  validation  of  a  kinetic 
scheme,  and  that  additional  comparisons  are  needed. 

This  then  leads  to  the  third  objective,  which  is  to  further  test  the  proposed  mechan¬ 
ism  against  a  wide  variety  of  experimental  data  found  in  the  literature  on  the  laminar 
flame  structure  obtained  from  flat-flame  burners,  and  on  the  homogeneous  reaction 
processes  obtained  from  flow  reactors,  static  reactors,  and  shock  tubes.  These 
additional  comparisons  cover  temperatures  from  820  to  2180K,  pressures  from  0.05 
to  4.7  atmospheres,  and  stoichiometries  from  very  fuel  lean  (<f>  =  0.05)  to  fuel 
pyrolysis  conditions.  It  will  be  shown  that  the  important  reactions  can  differ  for  the 
different  experimental  techniques,  thereby  demonstrating  the  importance  of  utilizing 
all  available  information  for  the  comprehensive  validation  of  a  proposed  kinetic 
mechanism,  as  originally  emphasized  by  Westbrook  and  Dryer  (1979). 

In  the  following  we  shall  first  define  the  experimental  and  numerical  aspects  of  the 
present  investigation.  This  will  be  followed  by  comparison  and  discussion  of  the 
results  from  the  various  studies  involving  the  laminar  flame  speeds,  the  structure  of 
burner-stabilized  flames,  the  flow  reactor,  the  static  reactor,  and  the  shock  tube. 


2  EXPERIMENTAL  METHODOLOGY 

The  counterflow  twin  flame  technique  for  the  determination  of  laminar  flame  speeds 
is  well  documented  (see  e.g..  Egolfopoulos  et  al .,  1989;  Zhu  et  al.,  1989).  It  involves 
the  establishment  of  two  symmetrical,  planar,  nearly-adiabatic  flames  in  a  nozzle¬ 
generated  counterflow'  configuration,  and  the  subsequent  determination  of  the  axial 
velocity  profile  along  the  centerline  of  the  flow  by:  laser  Doppler  velocimetry.  The 
minimum  point  of  the  velocity  profile  is  identified  as  a  reference  upstream  flame  speed, 
S„.  while  the  velocity  gradient  ahead  of  this  point  characterizes  the  imposed  stretch 
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rate.  K.  Thus  by  plotting  Su  vs  K,  the  stretch-free  flame  speed  S?  is  obtained  through 
linear  extrapolation  to  zero  stretch. 

The  main  difference  between  the  present  and  the  previous  investigation  is  that  a 
liquid  rather  than  a  gaseous  fuel  was  used,  hence  necessitating  the  inclusion  of  a 
continuous-flow  evaporator  which  vaporizes  and  mixes  the  methanol  with  preheated 
air.  The  air  was  filtered  and  preheated  by  passing  through  an  electrically-heated 
copper  tube,  while  methanol  vaporization  was  achieved  by  injecting  and  impinging  a 
liquid  jet  on  a  hot  plate.  All  flow  lines  were  electrically-heated  tn  order  to  avoid 
subsequent  condensation.  In  addition,  the  inner  surface  of  a  burner  was  heated  with 
a  continuous  flow  of  hot  water  through  a  surrounding  jacket.  The  unbumed  mixture 
temperature.  T„.  was  measured  with  a  thermocouple  at  the  exit  of  the  nozzle.  This 
temperature  was  controlled  by  the  heating  rates  of  the  air  and  the  internal  surfaces 
of  the  burners,  and  was  sufficiently  high  to  ensure  that  methanol  remains  in  the 
gaseous  phase  for  all  concentrations  within  the  flammable  range.  The  flow  rates  of 
methanol  and  air  were  respectively  measured  by  using  calibrated  flowmeters  and  sonic 
orifices  from  which  the  mixture  stoichiometry  could  be  determined.  The  stoichiometry 
was  further  verified  through  independent  gas  chromatographic  analyses  of  samples 
taken  at  the  exit  of  the  nozzle.  Agreements  were  obtained  for  all  cases.  Further  noting 
that  steady  state  was  maintained  during  the  several  hours  of  continuous  experimen¬ 
tation.  it  was  ensured  that  there  was  no  liquid  accumulation  in  the  flow  system. 


3  NUMERICAL  METHODOLOGY 

Numerical  simulation  of  all  systems  were  conducted  by  using  CHEMKIN-based 
programs  (Kee  et  a/..  1983;  Kee  et  al. ,  1985;  Lutz  et  al.,  1987;  Kee  et  al. ,  1989a)  and 
the  thermodynamic  data  base  developed  by  Kee  et  al.  (1987). 

For  the  laminar  flame  speed  determination,  the  freely-propagating  flame  version  of 
the  one-dimensional  code  (Kee  et  al..  1985;  Grcar  et  al.,  1986)  was  used,  while  for  the 
laminar  flame  structure  the  burner-stabilized  option  of  the  code  was  used.  Thermal 
diffusion  was  included  in  both  cases.  The  code  uses  a  hybrid  time-integration/ 
Newton-integration  technique  to  solve  the  steady-state  conservation  equations  of 
species  and  energy,  and  outputs  the  complete  flame  structure  and  the  sensitivities  of 
all  reaction  rates  on  the  dependent  parameters. 

For  the  simulation  of  flow  reactors,  static  reactors,  and  shock  tubes,  the  SENKIN 
code  (Lutz  et  al.,  1987)  was  used  which  predicts  homogeneous  gas  phase  chemical 
kinetics  by  integrating  the  time  dependent  conservation  equations  of  species  and 
energy.  The  code  was  also  modified  in  order  to  include  the  option  of  constant  volume 
and  temperature  reactions  which  is  relevant  to  static  reactor  studies. 

Detailed  information  on  the  use  of  the  codes  as  well  as  comments  on  the  validity 
of  the  comparisons  for  each  case  will  be  given  in  the  following  sections. 


4  STUDIES  BASED  ON  LAMINAR  FLAME  SPEEDS 

The  experimentally  determined  laminar  flame  speeds,  S°.  of  methanol/air  mixtures  at 
atmospheric  pressure  and  for  unburned  mixture  temperatures,  Tu,  of  318,  340,  and 
368  K  are  shown  in  Figure  2  for  a  wide  range  of  stoichiometry.  Figure  3  shows  the 
anticipated  result  that  S?  basically  increases  linearly  with  Ty  for  all  concentrations. 
Flame  speeds  at  298  K  were  not  measured  but  were  obtained  by  extrapolating  the  data 
at  higher  Tu  to  298  K.  As  such  they  should  be  treated  with  some  caution. 
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FIGURE  2  Experimentally  and  numerically  determined  laminar  flame  speeds.  S°(<t>.  TJ.  at  p  =  I  atm. 
using  the  present  kinetic  scheme. 

The  maximum  5"  for  methanol/air  is  found  to  occur  around  0  =  1.15,  which  is 
higher  than  the  corresponding  value  of  <t>  =  1.05  for  methane/air.  The  reason  for 
such  a  difference  is  the  existence  of  the  extra  O  atom  in  the  methanol  molecule  which 
provides  an  additional  amount  of  "oxidizer"  as  the  fuel  concentration  becomes  richer. 
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FIGURE  3  Experimentally  determined  laminar  flame  speeds.  Sj.  with  T„  for  lean,  stoichiometric,  and 
rich  methanol/air  flames  at  p  -  1  atm.  and  the  determination  of  S"  at  T„  =  298  K  through  linear 
extrapolation 
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FIGURE  4  Species  consumption  path  anaKsis  for  lean  methanol/air  flame  (<£  =  0.6)atT^  =  318Kand 
p  =  1  atm 

Furthermore  the  Sj's  for  methanol/air  mixtures  at  7"„  =  298  K  are  slightly  higher 
than  those  of  methane/air.  and  in  general  are  comparable  to  other  paraffin/air 
mixtures  at  the  same  initial  conditions. 

Figure  2  shows  the  numerically  calculated  results  and  their  close  agreements  with 
the  experimental  data.  The  kinetic  scheme,  listed  in  Table  I.  includes  30  species  and 
171  reversible  reactions  for  the  C, .  C:.  and  CH,OH  submechanisms.  Its  compilation 
was  conducted  iteratively  since  the  validity  of  certain  reactions  could  not  be  assessed 
b\  using  the  flame  speed  data  alone.  For  the  C,  and  C;  submechanisms,  the  starting 
point  was  the  mechanism  compiled  by  Egolfopoulos  et  al.  (1990.  1991)  (reactions  R1 
through  R 1 5 1 )  which  closely  predicts  the  laminar  flame  speeds  of  mixtures  of  CH4  and 
the  C: -hydrocarbons  with  0;  and  N:  for  various  stoichiometries  and  pressures.  For 
the  CH,OH  submechanism  a  combination  of  the  latest  information  of  Norton  and 
Dryer  (1989.  1990)  was  used  on  methanol  oxidation  and  pyrolysis  (reactions  R152 
through  R 1 7 1 ).  For  each  reaction  in  Table  I  the  forward  rate  is  reported  along  with 
its  source;  the  code  determines  the  backward  rate  through  thermodynamic  equilib¬ 
rium  using  the  thermnodynamic  data  of  Kee  et  al.  (1987).  The  fitted  backward  rates 
are  also  reported  in  Table  I  so  that  the  accurate  use  of  the  present  scheme  by  others 
is  facilitated. 

In  order  to  utilize  the  laminar  flame  speed  comparisons  for  the  validation  of  the 
kinetic  scheme,  it  is  important  to  identify  the  main  species  consumption  paths  in 
methanol  flames  and  conduct  sensitivity  analysis  (Kee  et  al.,  1985;  Grcar  et  al.,  1986) 
by  determining  the  influence  of  all  reaction  rates  on  the  laminar  mass  burning  rate 
m  =  o.,S[ '.  where  p_  is  the  density  of  the  unbumed  mixture. 

The  species  consumption  paths  are  shown  in  Figures  4.  5.  and  6  for  lean,  stoichio¬ 
metric.  and  rich  flames  respectively.  These  paths  were  determined  by  integrating  all 
reactions  through  the  flame  and  determining  the  fraction  of  each  species  consumed  by 
a  specific  reaction;  this  fraction  is  indicated  next  to  each  species.  From  such  a  study 
the  following  observations  can  be  made: 
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FIGURE  5  Species  consumption  path  analysis  for  stoichiometric  methanol /air  flame  (0  =  10)  at 
r  =  31KK  and  />  =  I  atm 

( I )  CH,OH:  It  is  attacked  by  H.  O.  OH.  and  HO;.  About  70  to  75%  is  converted  to 
the  hydroxymethyl  radical  (CH:OH).  while  only  25  to  30%  is  converted  to  the 
methoxy  radical  (CH,0).  Even  at  this  stage  the  very  different  nature  of  methanol 
chemistry  is  obvious  as  compared  to  methane  for  which  the  methyl  radical  (CH,)  is 
produced  after  the  initial  fuel  reactions.  For  mechanol  flames  it  is  clear  that  under¬ 
standing  of  the  kinetic  behavior  of  the  oxygenated  compounds.  CH;OH  and  CH,0. 


FIGURE  6  Species  consumption  paih  analysis  for  rich  methanol/air  flame  (0  =  20lai7'„  =  3l8Kand 
p  =  I  aim 
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is  crucial  to  the  compilation  of  a  reliable  kinetic  scheme,  while  the  knowledge  of  the 
CH,  chemistry  is  not  as  important.  Among  the  various  radicals  attacking  CH,OH,  the 
OH  consumes  almost  50-60%  for  lean  and  stoichiometric  flames  for  which  OH  is 
available,  and  decreases  to  about  25%  for  rich  flames  due  to  reduced  concentrations 
of  OH.  The  attack  by  H  is  less  important  for  lean  and  stoichiometric  flames  than  for 
rich  flames  where  it  contributes  to  almost  70%  of  the  fuel  consumption.  The  reaction 
with  O  also  becomes  progressively  less  important  as  the  mixture  becomes  richer  due 
to  the  decreased  O  concentration. 

(2)  CH,OH:  This  is  the  main  product  of  the  initial  fuel  reactions,  and  for  all  stoichi¬ 
ometries  there  are  only  two  consumption  paths.  The  most  important  is  the  reaction 
with  O,.  with  the  thermal  decomposition  route  becoming  somewhat  important  for 
rich  flames  due  to  the  reduced  supply  of  O, .  The  reaction  of  CH,OH  with  O,  produces 
CH:0  as  well  as  HO,.  Generation  of  the  latter  indicates  an  additional  difference 
between  methane  and  methanol  in  that  H02  can  also  be  produced  in  substantial 
amounts  by  means  other  than  reactions  R9  and  RIO  (H  +  O,  +  M  =  HO,  +  M). 
Therefore,  reactions  of  HO,  are  expected  to  be  of  increased  importance  as  compared 
to  those  in  hydrocarbon  flames. 

(3)  CH,0:  Its  main  consumption  paths  are  very  different  from  those  of  CH,OH, 
indicating  the  importance  of  considering  both  CH;OH  and  CH,0  in  recent  studies. 
The  methoxy  radicals  either  react  with  H  or  thermally  decompose.  Thermal  decom¬ 
position  leads  to  CH,0  and  it  is  more  important  for  lean  and  rich  flames  as  compared 
to  the  stoichiometric  flame.  The  reaction  with  H  is  the  main  source  of  CH,  radicals 
and  is  more  important  for  stoichiometric  flames  due  to  the  abundance  in  H.  For 
significantly  off-stoichiometric  mixtures,  the  H  concentration  decreases,  resulting  in 
an  increased  importance  of  the  thermal  decomposition  route.  It  is  interesting  to  note 
that,  in  contrast  to  CH,OH,  the  reaction  with  O,  does  not  appear  to  be  important. 

(4)  CH,:  It  exists  in  small  concentrations  because  it  is  produced  by  the  partial 
consumption  of  CH,0  whose  concentration  is  significantly  lower  than  that  of 
CH; OH.  For  lean  and  stoichiometric  flames  it  partially  reacts  with  O  leading  to 
CH,0;  for  all  stoichiometries  30-40%  of  the  CH,  reacts  with  HO,  leading  to  CH,0. 
The  reaction  with  HO;  is  important  due  to  the  increased  amount  of  H02  in  methanol 
flames.  The  reactions  with  H  and  CH,.  which  respectively  lead  to  CH4  and  C2- 
hydrocarbons.  become  more  important  for  rich  flames,  although  the  rates  of  these 
consumption  paths  are  lower  than  those  in  methane  flames. 

(5)  CH;0:  From  the  species  consumption  diagrams  shown  in  Figures  4,  5,  and  6,  it 
is  seen  that  for  all  stoichiometries  approximately  one  mole  of  CH,0  is  produced  from 
one  mole  of  CH,OH.  In  contrast,  in  methane  flames  the  direct  production  of  CH2  and 
HCO  radicals  as  well  as  C2  species  provides  alternate  paths,  as  shown  in  Figure  7 
where  the  species  consumption  paths  for  atmospheric,  stoichiometric  CH4/air  flames 
are  illustrated. 

A  comparison  between  methanol  and  methane  combustion  with  regard  to  the 
production  of  selected  species  is  shown  in  Figure  8.  The  reported  data  represent  the 
ratio  between  the  total  number  of  moles  of  each  species  present  throughout  a 
methanol/air  flame  and  the  corresponding  number  of  moles  for  a  methane/air  flame 
at  p  =  1  atm  and  <p  =  1  conditions.  Results  show  that  the  CH20  and  H02  concen¬ 
trations  present  in  a  methanol  flame  are  respectively  15  and  4  times  higher  than  those 
in  the  methane  flame.  The  concentrations  of  the  C2 -hydrocarbons  present  in  methanol 
flames  are.  however,  only  fractions  of  those  of  the  methane  flame.  These  results 
quantitatively  demonstrate  the  lower  propensity  of  methanol  to  soot,  due  to  the  lower 
C, -hydrocarbon  production,  and  at  the  same  time  its  higher  propensity  for  CHjO 
production. 
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FIGURE  8  Raticv  of  the  total  number  of  moles  of  selected  species  present  in  a  stoichiometric  methanol/ 
air  flame  (<f>  =  1.0)  at  Tu  -  298  K  and  p  =  1  atm  over  the  corresponding  number  for  stoichiometric 
methane/air  (</>  =  1.0)  flame  at  the  same  conditions.  Top:  CH:0.  HCO.  and  HO:.  Bottom:  C:-hydrocarbons 
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The  first  order  normalized  sensitivity  coefficients  (Kee  et  al .,  1985;  Grcar  et  al ., 
1986)  of  selected  reactions  on  the  mass  burning  rate.  m°,  of  lean,  stoichiometric,  and 
rich  flames  are  shown  in  Figure  9.  A  “shift”  of  the  importance  of  certain  reactions  can 
be  seen  as  the  stoichiometry  changes  from  lean  to  rich.  Based  on  the  results  of  Figures 
4.  5.  6  and  9.  the  role  of  the  key  reactions  on  the  overall  rate  progress  can  now  be 
assessed: 

(1)  Reaction  R1  (H  +  0:  =  OH  +  O):  It  has  high  positive  sensitivity  on  m°  for  all 
flames  due  to  its  chain  branching  nature,  and  is  the  most  dominant  reaction  for  rich 
flames.  Similar  behavior  exists  for  hydrocarbons  (Egolfopoulos  et  al.,  1990,  1991). 

(2)  Reactions  R9  (H  +  O,  +  M  =  HO,  +  M,  M  ±  H.O)  and  R10  (H  +  O,  + 
H,0  =  HO,  +  H,0):  They  compete  directly  with  R1  for  the  same  radicals,  pro¬ 
ducing  the  relatively  inactive  HO,  and  hence  providing  significant  termination  paths 
for  lean  flames.  For  rich  flames,  however,  the  termination  is  minimal. 

(3)  Reaction  R60  (CO  +  OH  =  CO,  -I-  H):  This  is  the  main  CO  oxidation  reaction, 
producing  heat  and  H  radicals.  For  lean  flames  its  positive  sensitivity  on  m°  dominates 
even  that  of  R1  because  there  is  enough  O,  to  oxidize  CO  to  CO,.  Its  importance, 
however,  is  gradually  diminished  as  the  stoichiometry  becomes  richer  because  there 
is  not  sufficient  O,  to  oxidize  CO.  As  a  result.  CO  assumes  the  role  of  an  abundant 
“final  product"  species  and  its  rate  of  oxidation  is  no  longer  a  controlling  factor. 

(4)  Reactions  R53  (HCO  +  M  =  CO  +  H  -I-  M)  and  R57  (HCO  +  Oj  = 
CO  +  HO;):  The  thermal  decomposition  reaction  R53  has  high  positive  sensitivity 
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FIGURE  10  Comparison  between  the  experimentally  determined  S?(d>)  of  methanol/air  flame  at 
Tu  =  368  K  and  numerical  calculations  using  the  present  scheme  with  the  rate  for  reaction  CH3OH  +  OH 
taken  from  Hess  and  Tully  (1989)  and  Bowman  (1975). 

on  m°  for  lean  flames  and  provides  termination  for  rich  flames.  Reaction  R57  is  an 
alternate  path  for  the  consumption  of  the  formyl  radicals  and  its  sensitivity  behaves 
exactly  opposite  to  that  of  R53.  For  lean  flames  its  rate  is  favored  over  R53  due  to 
the  excess  02  present,  and  provides  significant  termination  because  it  produces  H02 
which  is  less  essential  for  the  propagation  of  lean  flames  as  compared  to  H.  For  rich 
flames,  which  are  sensitive  to  all  potential  sources  of  OH,  R57  favors  propagation  by 
producing  additional  HO:  which  can  either  branch  to  20H  by  reacting  with  H,  or 
attack  the  fuel  molecule  through  R160  (CH3OH  +  H02  =  CH:OH  +  H202) 
which  produces  H20,  and  eventually  20H  through  decomposition.  As  a  result  its 
direct  competitor  R53,  which  produces  the  less  essential  CO  and  H  on  the  rich  side, 
has  a  negative  effect  on  m°.  For  laminar  flame  comparisons  the  rate  of  R57  was 
assumed  to  be  twice  that  of  Timonen  et  al.  ( 1 988),  in  which  the  uncertainty  factor  is 
stated  as  5.  This  modification  was  needed  in  order  to  obtain  closer  agreements  for  the 
laminar  flame  speeds  of  lean  mixtures,  which  were  initially  overpredicted  by  2-3  cm/s. 
(5)  Reactions  R157  (CH3OH  +  OH  =  CH2OH  +  H20)  and  R158  (CH3OH  + 
OH  =  CH,0  +  H20):  For  lean  flames,  propagation  is  favored  by  both  reactions 
but  the  associated  sensitivity  is  small  because  their  progress  is  fast,  leading  to  an 
almost  60%  fuel  consumption.  For  rich  flames,  however,  this  amount  drops  to  about 
25%  and  m°  becomes  more  sensitive  to  these  rates,  with  the  CH2OH  showing  a 
positive  sensitivity  while  the  competing  CH30  path  has  a  negative  effect.  The  sign 
reversal  of  the  sensitivity  of  the  CH30  path  for  stoichiometric  and  rich  concentrations 
is  due  to  the  fact  that  through  this  path  CH3,  CH4,  and  C2 -hydrocarbons  can  be 
formed  instead  of  CH;0.  This  can  retard  propagation  as  compared  to  the  CH2OH 
path  which  leads  directly  to  the  more  reactive  CH20.  The  rates  of  these  reactions  were 
obtained  from  two  sources.  The  first  was  from  Hess  and  Tully  (1989)  which  was  also 
used  by  Norton  and  Dryer  (1990).  The  second,  included  in  Table  I,  was  from  Bowman 
(1975)  with  a  rate  which  is  lower  than  that  of  Hess  and  Tully  by  a  factor  of  3.2  to  3.5 
for  T  >  1000  K.  Figure  10  shows  a  small  difference  in  the  calculated  flame  speeds 
using  these  two  rates.  Comparisons  with  the  turbulent  flow  reactor  oxidation  data  to 
be  discussed  subsequently,  however,  demonstrate  that  significant  differences  exist  in 
the  calculated  species  profiles  and  a  rate  which  follows  that  of  Bowman  rate  yields 
substantially  closer  agreement. 
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(6)  Reactions  R154  (CH,OH  +  H  =  CH2OH  +  H2)and  R155  (CH3OH  +  H  = 
CH,0  4-  H;):  They  are  important  for  rich  flames  with  theCH;OH  path  favoring  and 
the  CHjO  path  retarding  propagation,  for  reasons  similar  to  R157  and  R158 
respectively. 

(7)  Reaction  R160  (CH3OH  +  HO,  =  CH,OH  +  H,02):  It  appears  to  have  a 
significant  positive  effect  on  the  propagation  of  rich  mixtures.  Although  fuel  con¬ 
sumption  by  this  reaction  is  small,  approximately  4%,  its  sensitivity  is  high  because 
it  simultaneously  produces  two  species  whose  subsequent  reactions  favor  propa¬ 
gation.  That  is,  most  CH,OH  reactions  have  positive  sensitivity  on  m°  while  H202  can 
also  branch  to  20H  which  is  crucial  for  rich  flames.  This  demonstrates  that  both 
sensitivity  and  species  consumption  path  analyses  are  needed  for  accurate  assessment 
of  a  mechanism. 

(8)  Reactions  R 1 63  (CH2OH  +  M  =  CH20  -I-  H  +  M)and  R43  (CHjO  +  M  = 
CH,0  +  H  +  M):  Reaction  R163  does  not  have  a  significant  effect  on  propagation 
due  to  the  increased  CH,OH  consumption  through  R165  (CH2OH  +  02  = 
CH,0  +  HO,).  The  reported  rate  is  an  appropriately  fitted  one  reported  by  Tsang 
(1987).  Reaction  R43  has  a  positive  effect  on  propagation  for  all  stoichiometries  but 
its  importance  decreases  for  richer  flames.  The  rate  used  for  R43  is  the  one  reported 
by  Tsang  and  Hampson  (1986).  These  two  rates  are  not  well  established  and  their 
uncertainties  have  been  discussed  by  Grotheer  and  Kelm  (1989). 

(9)  HO,  Reactions.  Due  to  the  increased  production  of  HO,,  a  number  of  reactions 
are  shown  to  be  important  for  both  lean  and  rich  flames.  Reaction  R12 
(HO,  +  H  =  OH  4-  OH)  leading  to  20H  significantly  favors  propagation  for  all 
concentrations  and  its  direct  competitor  R 11  (HO,  4-  H  =  H,  4-  02),  leading  to  the 
stable  species  H2  and  02,  provides  significant  termination.  In  particular,  for  rich 
flames  Rll  is  the  dominant  termination  path.  Furthermore,  reactions  RI4 
(HO,  +  OH  =  H,0  4-  O,)  and  R15  (HO,  +  HO,  =  H,0,  +  O,)  have  a  notice¬ 
able  effect  on  propagation  as  well  as  the  H,0,  decomposition  reaction  RI6 
(H,0,  +  M  =  OH  +  OH  4-  M). 

During  compilation  of  the  present  scheme,  calculations  were  also  conducted  to 
assure  that  the  C,  and  C;  submechanisms  (reactions  R1  through  R151)  predict  the 
previously  reported  flame  speed  data  of  mixtures  of  CH4  and  the  C2 -hydrocarbons 
with  air  (Egolfopoulos  et  al .,  1989;  Egolfopoulos  el  al .,  1990,  1991).  Since  the 
thermodynamic  data  of  Kee  ei  al.  ( 1 987)  were  used  instead  of  the  data  of  Gordon  and 
McBride  (1971),  it  was  necessary  to  modify  the  forward  rate  of  reaction  R35 
(CH,  +  CH,  =  C,H5  +  H)which  was  shown  (Egolfopoulos  era/.,  1990, 1991)  to  be 
crucial  for  ethane  flames.  Specifically,  it  was  shown  that  the  backward  rate  of  R35  is 
an  important  termination  path  for  ethane  flames  and  that  in  order  to  closely  predict 
the  experimentally  determined  ethane/air  flame  speeds  a  much  faster  rate  was  needed 
than  those  reported  in  the  literature.  Moreover,  it  was  found  that  the  two  main 
sources  of  thermodynamic  data  (Gordon  and  McBride,  1971;  Kee  et  al.,  1987)  result 
in  very  different  equilibrium  constants  for  R35.  Values  obtained  from  Kee  ei  al.  (1987) 
are  10  times  higher  than  values  obtained  by  using  the  data  of  Gordon  and  McBride 
(1971)  in  the  temperature  range  of  1000  and  2500  K.  In  previous  modeling  (Egolfopoulos 
et  al..  1990.  1991)  the  forward  rate  of  Roth  and  Just  (1985)  and  the  low  equilibrium 
constant  of  Gordon  and  McBride  (1971)  were  used,  so  that  the  backward  rate  was 
sufficiently  hign  to  assure  agreements  with  the  experimental  flame  speed  data.  In  the 
present  investigation,  by  adopting  the  widely  accepted  data  of  Kee  et  al.  (1987),  a 
forward  rate  10  times  higher  was  used  for  R35.  In  order  to  achieve  closer  agreement  for 
the  laminar  flame  speeds  of  C,H4 /mixtures,  the  rate  of  Kee  et  al.  (1989b)  for  the  reaction 
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TABLE  111 

Comparison  between  experimental  hydrocarbon/air  laminar  flame  speeds  (cm/s),  Sj,  at  p  =  1  atm  and 
<p  =  1  with  predicted  values  using  the  present  kinetic  scheme 


CH, 

C,H6 

c,h4 

C2H, 

Experimental  data* 

40.2 

42.5 

65.8 

135.4 

Present  scheme 

39.3 

42.9b 

64.0b 

136.  O'1 

'As  determined  by  Egolfopoulos  el  al.  (1990.  1991). 
'’Calculated  with  C, -reactions  of  Table  11. 


RI06  (C:H,  +  H  =  C;H:  +  H,)  was  used  while  the  rate  of  reaction  R 1 10 
(C:H,  +  0:  =  CH20  +  HCO)  was  reduced  by  a  factor  of  four  (Egolfopoulos  et  al.. 
1990,  1991)  from  that  reported  by  Slagle  et  al.  (1989).  Furthermore,  the  formulation  of 
Stewart  et  al.  (1989)  for  the  pressure  dependence  (SRI  fall-off  parameters)  of  the 
reaction  R21  (CH4  +  M  =  CH,  +  H  +  M)  was  adopted.  Comparisons  between 
experimental  and  calculated  flame  speeds  for  stoichiometric  mixtures  of  CH4  and  the 
C2 -hydrocarbons  with  air  at  1  atm  are  shown  in  Table  III  and  close  agreement  can  be 
seen.  The  results  of  the  C;-hydrocarbons  were  obtained  by  also  including  the  C3- 
reactions  (R172  through  R219)  reported  in  Table  II.  C,  chemistry  has  been  shown 
(Egolfopoulos  et  al..  1990,  1991)  to  decrease  the  laminar  flame  speed  of  atmospheric, 
stoichiometric  C: -hydrocarbon/air  mixtures  by  1.5  to  2.5  cm/s. 

Grotheer  and  Kelm  (1989)  have  recently  conducted  a  detailed  numerical  study  of 
methanol  flames  with  the  latest  information  of  Norton  and  Dryer  (1989)  and  Tsang 
(1987).  predicting  laminar  flame  speeds  of  methanol/air  mixtures  at  p  =  1  atm  and 
Tu  =  298  K  conditions.  In  this  study,  the  one-dimensional,  time-dependent  flame  code 
of  Warnatz  (1978)  was  used  and  a  number  of  revisions  was  adopted  due  to  significant 
overprediction  of  the  flame  speeds  of  rich  mixtures.  The  revised  scheme  was  tested 
against  our  experimental  data  for  stoichiometric  flames  at  various  initial  temperatures 
using  the  Sandia  flame  code  (Kee  et  al.,  1985;  Grcar  et  al.,  1986),  and  the  results  are 
shown  in  Figure  1 1 .  First,  note  that  while  the  agreement  between  experimental  and 
numerical  data  is  good  at  low  temperatures,  the  temperature  dependence  is  not  captured 


FIGURE  II  Comparison  between  the  experimentally  determined  Sf's  of  stoichiometric  methanoi/air 
flame  (</>  =  l.Oiat  various  T,s  and  numerical  calculations  using  the  schemes  of  Grotheer  and  Kelm  (1989). 
Norton  and  Dryer  (1989.  1990).  the  C,  part  of  the  present  scheme,  and  the  “full"  present  scheme. 
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FIGURE  12  Comparison  between  the  experimental  S?'s  of  methanol/air  mixtures  as  determined  by 
Bradley  et  al.  (1991)  at  0.089  and  0.25  atm  and  numerical  calculations  using  the  present  scheme. 


closely,  leading  to  an  underprediction  of  about  6cm/s  for  the  Tu  =  368  K  flame. 
Second,  the  predicted  flame  speed  for  Tu  -  298  K  is  close  to  that  reported  by  Grotheer 
and  Kelm,  suggesting  that  the  numerical  codes  developed  by  Wamatz  and  Kee  and 
coworkers  lead  to  similar  results  for  the  same  kinetic  scheme.  It  may  also  be  noted  at 
this  point  that  during  the  compilation  of  the  present  kinetic  scheme  the  laminar  flame 
speeds  were  predicted  relatively  well  and  the  overpredictions  as  reported  by  Grotheer 
and  Kelm  were  not  encountered. 

In  Figure  1 1  the  present  experimental  data  are  also  compared  with  the  numerical 
predictions  by  using  the  two  schemes  reported  by  Norton  and  Dryer  (1989,  1990)  on 
methanol  oxidation  and  pyrolysis.  For  both  cases  significant  overpredictions  are  seen. 
Since  these  schemes  were  developed  based  mainly  on  flow  reactor  data,  the  observed 
discrepancy  demonstrates  the  importance  of  validating  proposed  schemes  against 
experimental  data  from  flames,  in  which  wide  variations  in  temperatures  and  concen¬ 
trations  are  embedded,  and  the  “coupling”  of  various  species  through  molecular 
diffusion  can  result  in  effects  which  cannot  be  adequately  captured  in  a  homogeneous 
system. 

In  order  to  quantify  the  contribution  of  the  C2  species  chemistry  on  laminar  flame 
propagation,  calculations  were  conducted  for  the  stoichiometric  methanol/air  flames  at 
various  Tu' s,  both  with  ?nd  without  the  C2  chemistry  of  the  present  scheme.  These 
results  are  compared  with  the  experimental  data  in  Figure  1 1 ,  where  it  can  be  seen  that 
the  effect  of  the  C2  chemistry  is  to  slightly  reduce  the  flame  speed,  by  3-4  cm/s.  This  effect 
is  due  to  the  fact  that  with  the  inclusion  of  the  C2  chemistry,  the  CH3  radicals  partially 
recombine  to  form  C2H6  instead  of  being  oxidized  to  CH20  or  CH30  which  are  less 
stable. 

The  present  scheme  was  further  tested  against  the  new  experimental  flame  speed 
data  as  determined  by  Bradley  et  al.  (1991)  at  0.089  and  0.25  atm.  The  comparisons 
are  shown  in  Figure  1 2  and  it  can  be  seen  that  although  the  calculations  predict  closely 
the  0.089  atm  data,  they  overpredict  the  flame  speeds  at  0.25  atm.  The  calculations 
performed  by  Bradley  et  al.  ( 1 99 1 )  predicted  closely  the  0.25  atm  data  but  significantly 
underpredicted  the  data  at  0.089  atm,  indicating  that  their  scheme  is  slower  than  the 
present  one.  Since  Bradley  et  al.  (1991)  mentioned  that  there  were  some  uncertainties 
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FIGURE  13  Comparison  between  the  experimental  (symbols)  flame  structure  as  determined  by 
Vandooren  and  Van  Tiggelen  (1981)  and  numerical  calculations  (lines)  using  the  present  scheme.  (Flame  1: 
m"  =  3 .302E  —  03 gm/sec-cnr .  <t>  =  0.89.  19.9%  CH,OH.  33.7%  0:.  46.4%  Ar.  p  =  0.053atm;  all  con¬ 
centrations  on  per  mole  basis). 

associated  with  the  measurements,  further  attempts  to  assess  the  cause  of  these  dis¬ 
crepancies  were  not  made. 

5  STUDIES  BASED  ON  LAMINAR  FLAME  STRUCTURES 

Since  the  laminar  flame  speed  is  only  a  global  property  of  the  flame,  its  sensitivity  on 
certain  aspects  of  the  flame  kinetics  can  be  quite  limited.  A  more  stringent  test  would 
involve  predictions  and  comparisons  of  the  flame  structure,  especially  the  stable  and 
radical  species  profiles.  Among  the  various  experimental  studies  which  have  been 
conducted  on  low  pressure  methanol  flames  using  the  flat-flame  burner  technique,  the 
experimental  data  reported  by  Vandooren  and  Van  Tiggelen  (1981),  Pauwels  et  al. 
(1989.  1990),  and  Bradley  et  al.  (1991)  were  modeled  herein. 

In  the  study  by  Vandooren  and  Van  Tiggelen  (1981),  the  molecular  beam  sampling 
technique,  coupled  with  mass  spectrometric  detection,  were  used  to  measure  the 
concentrations  of  both  stable  species  and  radicals  such  as  H,  O,  and  OH.  Conse¬ 
quently,  rate  coefficients  were  derived  for  the  reactions  of  (a)  methanol  with  H  and 
OH,  (b)  hydroxymethyl  with  oxygen,  and  (c)  the  thermal  decomposition  of  for¬ 
maldehyde.  In  the  studies  by  Pauwels  et  al.  (1989, 1990),  sampling  was  performed  with 
a  quartz  probe  and  the  concentrations  were  determined  by  using  electron  spin 
resonance  spectroscopy  for  H,  O,  and  OH,  and  gas  chromatography  for  the  stable 
species.  These  data  were  tested  (Pauwels  et  al.,  1990)  against  calculations  using  an 
updated  version  of  the  Westbrook  and  Dryer  mechanism  (1979)  with  34  reversible 
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FIGURE  14  Comparison  between  the  experimental  (symbols)  flame  structure  as  determined  by 
Vandooren  and  Van  Tiggelen  ( 1981 )  and  numerical  calculations  (lines)  using  the  present  scheme.  (Flame  II: 
m"  =  3.029  x  10  'gm/sec-cm;.  4>  =  0.36. 19.4%  CH, OH.  80.6%  0,.p  =  0.053  atm:  all  concentrations  on 
per  mole  basis). 

reactions,  and  it  was  found  that  the  structure  of  fuel  lean  flames  was  not  reproduced 
accurately.  Bradley  et  al.  (1991)  reported  the  structure  of  a  moderately-rich  flame 
determined  by  microprobe  for  species  sampling,  gas  chromatography  for  species 
analysis,  and  LDV  for  the  flow  velocities.  A  kinetic  scheme  was  also  compiled  which 
closely  reproduced  their  measured  flame  structure  but  underpredicted  the  corre¬ 
sponding  flame  speed  by  1 1  cm/s.  In  all  three  studies  (Vandooren  and  Van  Tiggelen, 
1981;  Pauwels  et  al.,  1989,  1990;  and  Bradley  et  al..  1991),  temperature  was  measured 
with  thermocouple  and  corrected  for  radiative  losses.  The  measurements  of  Pauwels 
et  al.  (1989.  1990)  had  an  additional  conductive  loss  to  the  sampling  probe  by  being 
adjacent  to  the  thermocouple. 

For  the  modeling  of  burner-stabilized  flames,  the  Sandia  code  (Kee  et  al.,  1985; 
Grcar  et  al..  1986)  was  used  which  can  either  solve  the  energy  equation  or  employ  a 
user-specified  temperature  profile.  Experimentally,  measurements  close  to  the  burner 
surface  can  be  difficult  to  obtain  and  data  are  typically  reported  only  beyond  a  certain 
minimum  distance  from  the  burner.  For  example,  in  the  studies  by  Vandooren  and 
Van  Tiggelen  (1981)  and  Pauwels  et  al.  (1989.  1990).  the  first  experimental  concen¬ 
trations  were  reported  at  temperatures  of  ~  1000  K.  Moreover,  the  use  of  a  thermo¬ 
couple  for  temperature  measurements  can  result  in  locally  reduced  temperatures  and 
that  imposes  the  constraint  of  using  the  experimental  temperature  profile  instead  of 
solving  the  energy  equation.  This  approach  both  facilitates  the  calculations  and  more 
accurately  simulates  the  actual  experiment  since  heat  losses  do  not  have  to  be 
modeled.  It  must  be  recognized,  however,  that  by  not  solving  the  energy  equation  the 
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FIGURE  15  Comparison  between  the  experimental  (symbols)  flame  structure  as  determined  by 
Vandooren  and  Van  Tiggelen  (1981)  and  numerical  calculations  (lines)  using  the  present  scheme.  (Flame  Ill: 
m°  =  5.207  *  10_,gm/sec-cnr .  <t>  =  0.21. 10.9%  CH, OH.  3.2%  H,.  85.9%  02,p  =  0.053 atm;  all  concen¬ 
trations  on  per  mole  basis). 

exothermicities  of  certain  reactions  are  imposed  instead  of  being  independently 
calculated.  Furthermore,  by  using  intrusive  techniques  for  flame  structure  studies 
there  is  always  the  possibility  of  altering  local  species  concentrations  due  to  inter¬ 
actions  with  the  probe. 

In  the  present  investigation,  the  Sandia  code  was  used  with  the  experimental 
temperature  profile  as  input  for  the  flames  of  Vandooren  and  Van  Tiggelen  (1981)  and 
Pauwels  el  al.  (1989,  1990).  The  data  of  Bradley  et  al.  (1991)  were  simulated  by  also 
solving  for  the  energy  equation  since  data  were  reported  close  to  the  burner  surface 
and  the  authors  have  also  modeled  the  flame  as  adiabatic.  For  all  cases,  at  the  first 
computational  point  the  summation  of  the  species  mole  fractions  ranges  between  0.95 
and  0.99.  indicating  that  all  species  were  mostly  accounted  for. 

Figures  13,  14..  and  15  show  the  comparisons  between  calculations  using  the 
present  scheme  and  the  experimental  data  for  flames  I,  II,  and  III  respectively,  as 
reported  by  Vandooren  and  Van  Tiggelen  (1981),  at  p  =  0.053  atm  (40Torr).  Flame  I 
is  a  lean  flame  (<p  =  0.89)  with  molar  concentrations  of  19.9%  CH3OH,  33.7%  02, 
and  46.4%  Ar,  and  the  comparisons  are  shown  in  Figure  13.  It  is  seen  that  in  the 
methanol  consumption  zone  the  profiles  of  all  stable  species  are  predicted  closely,  with 
the  exception  of  02  which  is  underpredicted.  The  profiles  of  H,  O,  and  OH  are 
predicted  satisfactorily,  but  due  to  uncertainties  associated  with  the  measurement  of 
these  radicals  it  is  difficult  to  rigorously  assess  the  validity  of  the  predictions.  Flame 
II  is  a  lean  flame  (<f>  =  0.36)  with  molar  concentrations  of  19.4%  CH3OH  and  80.6% 
02 .  and  the  comparisons  are  shown  in  Figure  1 4.  The  majority  of  stable  species  as  well 
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Spatial  Distance  X,  cm 

FIGURE  16  Comparison  between  the  experimental  (symbols)  flame  structure  as  determined  by  Pauwels 
et  al.  (1989.  1990)  and  numerical  calculations  (lines)  using  the  present  scheme.  (m°  =  2.96  x  I0~5gm/ 
sec-cm\  0  =  0.82.  9.9%  CH,OH,  18.1%  02.  72.0%  N,,  p  =  O.IOSatm;  all  concentrations  on  per  mole 
basis). 


as  H.  O.  and  OH  are  predicted  closely.  The  exceptions  are  H;,  which  exists  in  trace 
amounts,  and  CH,OH.  whose  predicted  rate  of  consumption  is  somewhat  faster  than 
the  experimental  values.  Flame  III  is  a  lean  flame  (<j>  =  0.21)  with  molar  concen¬ 
trations  of  10.9%  CHiOH,  3.2%  H:.  and  85.9%  0;,  and  the  comparisons  are  shown 
in  Figure  15.  All  stabie  species  are  predicted  closely,  but  the  H,  O,  and  OH  radicals 
are  significantly  underpredicted.  The  close  agreements  of  the  CHjOH  profile  for 
Flames  I  and  III  contrast  with  the  discrepancy  found  in  Flame  II,  suggesting  that 
systematic  errors  in  the  methanol  consumption  reactions  cannot  be  identified  through 
these  comparisons. 

Figures  16.  17,  and  18  compare  calculations  using  the  present  scheme  and  the 
experimental  data  as  reported  by  Pauwels  et  al.  (1989,  1990)  for  methanol/air  flames 
at  p  =  O.IOSatm  and  equivalence  ratios  of  0.82,  1.08,  and  1.50  respectively.  For  all 
cases  the  stable  species  are  predicted  closely  within  the  fuel  consumption  zone,  with 
the  exception  of  the  CHjOH  profile  for  the  4>  -  1.50  flame  which  is  “steeper”  as 
compared  to  that  of  the  experiment.  The  radical  species  H,  O,  and  OH,  however,  are 
overpredicted  and  this  disagreement  is  more  significant  for  the  H  radicals  in  the 
6  =  1.50  flame.  This  is  probably  the  cause  for  the  faster  consumption  rate  of  CH3OH 
through  its  reactions,  R154  and  R155.  with  H.  The  location  at  which  radicals  start 
being  produced  is  predicted  closer  as  compared  to  the  calculations  of  Pauwels  et  al. 
(1990). 

Figure  19  compares  calculations  using  the  present  scheme  and  the  experimen¬ 
tal  data  as  reported  by  Bradley  et  al.  (1991)  for  a  methanol/air  flame  at  p  = 
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FIGURE  1 7  Comparison  between  the  experimental  (symbols)  flame  structure  profiles  as  determined  by 
Pauwels  ei  a/  (1989  1990)  and  numerical  calculations  (lines)  using  the  present  scheme. 

<m°  =  3.15  *  lO-’gm'sec-cm-'.  <t>  =  1 .08.  12.5%  CH.OH.  17.5%  0:.  70.0%  N;.  p  =  0  1 05 atm:  all 
concentrations  on  per  mole  basis). 


0.089  atm  and  equivalence  ratio  4>  =  1.25.  Results  show  a  very  close  agreement  for 
all  species  with  the  exception  of  some  overprediction  of  the  CO  after  its  peak.  The 
present  scheme  also  reproduced  very  closely  both  the  temperature  and  flow  velocity 
profiles. 

For  the  comparisons  conducted  in  Figures  13  through  19,  a  general  conclusion 
is  that  the  present  scheme  appears  to  closely  predict  the  flame  structure  in  most 
cases,  capturing  the  most  significant  trends  of  the  profiles  of  the  stable  species 
concentrations.  In  particular,  the  predicted  CO  profiles  are  very  close  to  the  experi¬ 
mental  data  both  in  magnitude  and  also  in  the  location  of  the  maximum.  Recog¬ 
nizing  that  CO  is  an  important  intermediate  and  its  oxidation  can  significantly 
influence  global  flame  properties  such  as  the  laminar  flame  speed,  it  is  reasonable 
that  the  accurate  prediction  of  laminar  flame  speeds  by  the  present  kinetic  scheme 
for  various  conditions  correctly  accounts  for  the  CO  profile.  In  this  study,  attempts 
to  obtain  better  agreements  for  the  H,  O,  and  OH  radicals  were  not  made.  Experi¬ 
mental  measurement  of  these  radicals  is  much  more  difficult  than  that  of  the 
stable  species  because  they  exist  in  very  small  amounts  and  they  are  more  susceptible 
to  recombinations  when  intrusive  probing  techniques  are  used.  Further  improve¬ 
ment  in  studies  of  this  nature  would  require  non-intrusive  laser  diagnostic  techniques 
to  provide  concentration  and  temperature  results  with  higher  degrees  of  confidence 
and  accuracy. 

Finally,  calculations  of  the  flame  structures  studied  herein  were  made  by  using  the 
present  scheme  and  the  two  different  rates  (Bowman,  1975;  Hess  and  Tully,  1989)  for 
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FIGURE  18  Comparison  between  the  experimental  (symbols)  flame  structure  as  determined  by  Pauwels 
ei  at.  (1989.  1990)  and  numerical  calculations  (lines)  using  the  present  scheme.  (m°  =  2.79  x  10'Jgm/ 
sec-cnr.  <t>  =  1.50.  14.6%  CH,OH.  16.6%  O,.  66.8%  N;.  p  =  0.105atm;  all  concentrations  on  per  mole 
basis). 

R157  (CH,OH  +  OH  =  CH2OH  +  H;0)  and  R158  (CH,OH  +  OH  =  CHjO  + 
H:0)  reactions.  For  all  cases  it  was  found  that  the  CH,OH  profile  was  insensitive  to 
the  choice  of  the  rate  for  R157  and  R158,  indicating  an  additional  difference  between 
flame  studies  and  pure  kinetic  studies  as  it  will  be  further  demonstrated  in  the  next 
section. 


6  STUDIES  BASED  ON  FLOW  REACTOR  RESULTS 

The  turbulent  flow  reactor  is  a  valuable  experimental  technique  for  the  understanding 
and  development  of  kinetic  schemes.  Provided  that  there  are  no  steep  concentration 
gradients  in  the  system,  molecular  diffusion  is  negligibly  small  compared  to  axial 
convection  and  therefore  the  reaction  progress  is  solely  kinetically  controlled.  The 
axial  distance  can  be  directly  converted  to  the  time  coordinate  through  the  local  flow 
velocity,  .and  modeling  is  conducted  by  solving  the  time  dependent  equations  of 
species  and  energy.  The  flow  reactor  operates  in  temperatures  between  approximately 
1000  and  1200K,  and  the  system  is  near-adiabatic  before  any  significant  amount  of 
heat  is  released.  Therefore,  for  the  numerical  modeling  of  the  flow  reactor  data, 
adiabatic  solutions  were  adopted  so  that  the  exothermicity  and  endothermicity  of 
individual  reactions  are  also  accounted.  This  procedure  was  necessary  especially  for 
the  oxidation  studies  in  which  significant  heat  is  released.  In  the  pyrolysis  studies,  the 
system  can  be  considered  as  practically  isothermal  since  CO  is  not  oxidized.  In  the 
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FIGURE  19  Comparison  between  the  experimental  (symbols)  flame  structure  as  determined  by  Bradley 
ei  al.  (1991)  and  numerical  calculations  (lines)  using  the  present  scheme.  (m°  =  6.41  x  10'  3gm/sec-cm2, 
<t>  =  1.25.  14.9%  CH, OH.  17.88%  O..  67.22%  N,.p  =  0.089  atm;  all  concentrations  on  per  mole  basis). 


present  investigation  the  data  of  Norton  and  Dryer  (1989)  and  Aronowitz  et  al.  (1977, 
1979)  for  oxidation  and  pyrolysis  conditions  were  simulated. 

It  may  be  recalled  from  the  laminar  flame  speed  studies  that  modifications  of  the 
methanol  submechanism  were  not  required  since  these  reactions  do  not  significantly 
affect  propagation.  An  earlier  version  of  the  present  mechanism  was  tested  against  the 
oxidation  data  of  Norton  and  Dryer  (1989)  by  using  SENKIN  (Lutz  et  al .,  1987)  for 
a  constant  pressure,  adiabatic  system  which  closely  resembles  the  experimental  con¬ 
figuration.  In  these  simulations,  a  discrepancy  was  found  between  numerical  and 
experimental  results  as  indicated  by  a  predicted  methanol  profile  which  was  too 
“steep"  and  a  temperature  increase  which  was  substantially  faster  as  compared  to  the 
experimental  data.  Similar  behavior  was  also  observed  when  modeling  the  oxidation 
data  of  Aronowitz  et  al.  (1979). 

The  observed  discrepancies  motivated  further  investigation  of  the  methanol  sub¬ 
mechanisms  of  Norton  and  Dryer  (1989,  1990),  and  the  sensitivity  of  each  reaction 
of  CH,OH.  CH;OH.  CH,0.  CH;0,  HCO.  HO,,  and  H,0,  to  the  fuel  decay  profile 
was  examined  in  detail.  These  species  were  chosen  because  they  are  the  important 
intermediates  subsequent  to  fuel  decay,  and  because  they  can  also  significantly 
influence  the  production  of  H  and  OH  radicals.  It  was  consequently  identified  that  the 
methanol  profile  is  primarily  sensitive  to  the  reactions  between  CH3OH  and  OH 
which  produce  CH,OH  and  CH,0  (reactions  R157  and  R158  respectively).  By 
examining  the  reactions  controlling  the  OH  production,  it  was  found  that  by  increas¬ 
ing  the  rate  of  R1 1  (HO;  +  H  =  H2  +  02)  by  a  factor  of  3,  improved  agreements 
with  the  experimental  methanol  profile  were  obtained.  Reaction  Rll  provides  ter- 
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FIGURE  20  Comparison  between  the  experimental  (symbols)  profiles  for  CH,OH  at  various  <p’s  as 
determined  by  Norton  and  Dryer  (1989)  in  a  flow  reactor,  and  numerical  calculations  (lines)  using  the 
present  scheme  with  the  rate  for  reaction  CH,OH  +  OH  taken  from  Hess  and  Tully  (1989)  and  Bowman 
(19751. 


mination  to  the  system  since  it  does  not  allow  HO;  to  further  branch  to  OH.  This 
modification,  however,  was  not  adopted  because  certain  important  intermediates  such 
as  CH;0  and  H:  were  overpredicted  by  as  much  as  a  factor  of  two,  while  the  02  profile 
was  not  reproduced  satisfactorily  and  the  temperature  profile  also  appeared  to  have 
a  much  wider  “plateau”  (Norton  and  Dryer,  1989)  as  compared  to  the  experiments. 
Furthermore,  it  has  been  shown  (Egolfopoulos  and  Law,  1991)  that  the  low  tem¬ 
perature  hydrogen  chemistry  is  slow  in  predicting  the  laminar  flame  speeds  of  very 
lean  hydrogen/air  mixtures  and  this  modification  would  further  retard  the  overall 
reaction  progress.  The  rates  of  Hess  and  Tully  (1989)  for  reactions  R157  and  R158, 
used  by  Norton  and  Dryer  (1990),  were  measured  up  to  866  K  and  are  higher  than 
previously  recommended  rates  (Bowman,  1975;  Westbrook  and  Dryer,  1979)  for 
temperatures  above  700  K.  The  branching  ratio  CH:OH :  CH,0  used  by  Norton  and 
Dryer  (1990)  was  1  :  1,  as  suggested  by  Hess  and  Tully.  and  in  the  first  attempt  to 
improve  on  the  predictions  this  ratio  was  varied  from  0.2 ;  0.8  to  0.8 : 0.2  with  the  total 
rate  being  held  constant.  That  had  an  effect  on  the  induction  time  but  the  slope  of  the 
methanol  profile  remained  practically  unchanged.  The  reactor  temperature  was  also 
varied  by  ±  I5K,  resulting  in  only  minor  changes  in  the  predicted  slope  of  the 
methanol  profile.  Finally,  it  was  found  that  a  total  rate  for  R157  and  R158  smaller 
by  3.5  times  at  ~  1000  K  would  result  in  very  close  agreements  for  the  methanol 
profiles  as  well  as  the  other  species  and  temperature.  Noting  that  the  rate  suggested 
by  Bowman  (1975)  for  high  temperatures  appears  to  be  almost  3.2  times  slower  at 
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FIGURE  21  Comparison  between  the  experimental  (symbols)  profiles  for  CH3OH  at  various  <t>' s  as 
determined  by  Aronowiiz  ei  at.  ( 1979)  in  a  flow  reactor  and  numerical  calculations  (lines)  using  the  present 
scheme  with  the  rates  for  reaction  CH,OH  +  OH  taken  from  Hess  and  Tully  (1989)  and  from  Bowman 
(1975). 

1000  K.  Bowman's  rate  and  a  branching  ratio  of  1  : 1  were  subsequently  adopted  and 
extensive  agreements  were  obtained  for  all  oxidation  data. 

In  order  to  demonstrate  the  effect  of  these  modifications,  the  experimental  data  for 
methanol  decay  in  all  oxidation  runs  studied  herein  are  compared  with  the  predictions 
of  Figures  20  and  21  by  using  the  two  rates  for  R157  and  R158.  It  is  clear  that 
significant  improvement  for  lean  and  rich  cases  is  obtained  by  using  the  slower  rate. 
The  scheme  of  Grotheer  and  Kelm  ( 1 989)  was  also  tested  against  the  flow  reactor  data 
and  it  was  found  that  it  results  in  much  longer  induction  times,  as  compared  to  the 
present  scheme,  as  well  as  unsatisfactory  species  predictions. 

Predictions  from  calculations  by  using  Bowman's  rate  (1975)  for  R157  and  R158 
are  shown  in  Figures  22  and  26  for  various  experimental  data  sets  on  oxidation.  In 
all  cases,  due  to  the  uncertainty  related  to  the  experimental  induction  time,  the 
numerical  results  had  to  be  "time  shifted”  (Norton  and  Dryer,  1989,  1990;  Norton, 
1990)  in  order  to  match  the  experimental  data  at  the  50%  fuel  decay  point;  this 
amount  is  reported  in  each  figure.  The  results  show  that  there  is  good  agreement  for 
most  of  the  species  until  the  temperature  increases  to  the  point  where  heat  loss  to  the 
wall  (maintained  at  ~  1000  K)  and/or  molecular  diffusion  become  significant.  In  this 
regime  the  adiabatic,  homogeneous  model  is  probably  not  a  good  representation  of 
the  experimental  conditions  and  thus  discrepancies  in  the  CO  oxidation  region 
immediately  following  the  fuel  consumption  zone  were  not  further  assessed.  For  all 
cases,  however,  the  calculated  peak  values  of  CO  as  well  as  their  locations  agree  very 
well  with  the  experimental  data.  In  addition,  the  present  mechanism  captures  the 
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FIGURE  22  Comparison  between  experimental  (symbols)  flow  reactor  oxidation  data  for  <t>  -  0.6  as 
determined  bv  Norton  and  Dryer  (1989)  and  numerical  calculations  (lines)  using  the  present  scheme  of 
Table  I.  <0.779%  CH, OH.  1.98%  0:.  97.2404%  N,.  T,„  =  1027K,p,„  =  1. 0 atm,  adiabatic,  isobaric;  all 
concentrations  on  per  mole  basis). 


temperature  “plateau”  for  rich  flames  as  well  as  the  “early”  C02  production,  as  first 
analyzed  by  Norton  and  Dryer  (1989).  It  is  also  important  to  note  the  close  prediction 
of  the  concentration  profiles  of  CH:0,  which  is  a  key  species  for  methanol  oxidation. 

The  present  mechanism  has  been  also  tested  against  the  pyrolysis  flow  reactor  data 
of  Aronowitz  el  al.  (1977)  and  the  results  are  shown  in  Figure  27.  The  numerical 
solution  was  obtained  for  isothermal,  constant  pressure  conditions.  The  agreements 
for  all  species  are  very  close;  in  particular,  the  CH4  profile  is  predicted  closer  than  that 
of  Norton  and  Dryer  ( 1 990),  possibly  due  to  the  completeness  of  the  C,  and  C2  subsets 
of  the  present  mechanism. 


FIGURE  23  Comparison  between  experimental  (symbols)  flow  reactor  oxidation  data  for  <)>  =  1.2  as 
determined  by  Norton  and  Dryer  (1989)  and  numerical  calculations  (lines)  using  the  present  scheme  of 
Table  I.  (0.943%  CH.OH.  1.16%  0;.  97.8964%  N;.  Tm  =  1030 K.  />,„  =  I.Oatm.  adiabatic,  isobaric;  all 
concentrations  on  per  mole  basis). 
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FIGURE  24  Comparison  between  experimental  (symbols)  flow  reactor  oxidation  data  for  <t>  -  1.6  as 
determined  bv  Norton  and  Dryer  (1989)  and  numerical  calculations  (lines)  using  the  present  scheme  of 
Table  I.  (1.01%  CH.OH.  0.961%  O;.  98.0284%  N,.  T„  =  1034 K.  />,„  =  l.Oatm.  adiabatic,  isobanc:  all 
concentrations  on  per  mole  basis). 

7  STUDIES  BASED  ON  STATIC  REACTOR  RESULTS 

In  order  to  assess  the  validity  of  the  present  kinetic  scheme  at  lower  temperatures, 
static  reactor  data  were  also  modeled.  An  important  kinetic  difference  between  the 
static  reactor  and  other  experimental  configurations  is  that  lower  initial  mixture 
temperatures  can  be  obtained  since  there  is  no  geometric  restrictions  for  a  long 
induction  time.  Therefore,  the  characteristic  time  scale  for  the  static  reactor  (sec)  can 
be  orders  of  magnitude  greater  than  the  typical  time  scales  for  either  the  flow  reactor 
(msec)  or  the  shock  tube  Otsec).  Experimental  studies  of  methanol  pyrolysis  and 
oxidation  in  static  reactors  have  been  conducted  by  Cathonnet  er  al.  (1979.  1982)  for 
various  pressures  and  temperatures;  for  all  experiments  the  reactor  is  reported  as  near 


FIGURE  25  Comparison  between  experimental  (symbols)  flow  reactor  oxidation  data  for  <f>  =  0.05  as 
determined  by  Aronowitz  el  at.  (1979)  and  numerical  calculations  (lines)  using  the  present  scheme  of 
Table  I.  (0.6754%  CH,OH.  20.2613%  O..  79.0633%  N';.  7]„  =  1000K.pin  =  l.Oatm,  adiabatic,  isobaric: 
all  concentrations  on  per  mole  basis). 
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FIGURE  26  Comparison  between  experimental  (symbols)  flow  reactor  oxidation  data  for  <t>  =  1.70  as 
determined  by  Aronowitz  el  al.  (1979)  and  numerical  calculations  (lines)  using  the  present  scheme  of 
Table  1.  <0.6601%  CH,OH.  0.5824%  02,  98  7575%  NJt  Tm  «=  998K,p,„  =  I.Oatm,  adiabatic,  isobaric; 
all  concentrations  on  per  mole  basis). 


isothermal.  The  modeling  was  conducted  for  constant  volume  and  temperature 
conditions  in  order  to  closely  simulate  the  reported  experimental  configuration,  and 
in  no  case  was  it  necessary  to  “time  shift”  the  numerical  results  in  order  to  match  the 
experimental  data. 

Two  sets  of  the  pyrolysis  data  from  the  study  of  Cathonnet  et  al.  (1979)  were 
modeled  for  the  conditions:  (a)  p  =  0.263  atm  and  T  =  918  K  and  (b)p  =  0.526  atm 
and  T  =  973  K.  In  both  cases  the  inert  gas  was  nitrogen.  Comparisons  between 
calculated  and  experimental  data  are  shown  in  Figure  28  and  the  agreement  is  within 
the  reported  temperature  uncertainty  of  ±  1%.  Norton  and  Dryer  (1990)  conducted 
a  sensitivity  analysis  on  the  relative  effects  of  separately  increasing  the  rate  of  reaction 
R170  (CH;OH  +  CH2OH  =  CH3OH  +  CH20)  by  a  factor  of  5  and  varying  the 
temperature  by  10  K.  It  was  shown  that  both  adjustments  can  have  similar  effects  on 


FIGURE  27  Comparison  between  experimental  (symbols)  flow  reactor  pyrolysis  data  as  determined  by 
Aronowitz  et  al.  (1977)  and  numerical  calculations  (lines)  using  the  present  scheme  of  Table  I.  (0.815% 
CH,OH.  99.19%  N,.  Tm  =  I.Oatm.  isothermal,  isobaric;  all  concentrations  on  per  mole  basis). 
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FIGURE  28  Comparison  between  experimental  (symbols)  static  reactor  pyrolysis  data  as  determined  by 
Cathonnet  ei  at.  (1979)  and  numerical  calculations  (lines)  using  the  present  scheme  of  Table  I.  (Top:  1 1 .6% 
CH-.OH.  88.4%  N,.  Tm  =  918K.  pm  =  0.263 atm,  isothermal,  isochoric:  Bottom:  5.9%  CH,OH.  94.0% 
N;.  T,„  =  973  K.  p,„  =  0.526  atm.  isothermal,  isochoric;  all  concentrations  on  per  mole  basis). 

the  methanol  profile  and  therefore  rate  modifications  were  not  made.  The  present 
investigation  was  guided  by  this  approach  and  further  investigation  was  conducted. 
It  was  found  that  the  increase  of  reaction  R1 70  by  a  factor  of  5  provides  in  some  cases 
closer  agreements  for  the  methanol  profile.  In  general,  however,  better  agreements  for 
the  intermediates  such  as  H2  and  CO  were  provided  by  using  the  rate  which  is  reported 
in  Table  I  and  by  varying  the  temperature  within  its  uncertainty  range.  Since  tem¬ 
perature  variations  alone  can  account  for  improved  agreement  between  experimental 
and  numerical  results,  there  is  insufficient  justification  for  the  increase  of  the  rate  of 
this  reaction  from  the  value  reported  by  Norton  and  Dryer  (1990).  However,  further 
investigation  of  this  reaction  is  recommended  due  to  its  importance  under  fuel 
pyrolysis  conditions. 

In  Figures  29  and  30  experimental  oxidation  data  from  the  study  of  Cathonnet  ex 
al.  (1982)  are  compared  with  predictions  and  satisfactory  agreements  can  be  seen.  The 
system  temperature  was  varied  by  ±  1  %  in  order  to  investigate  .he  sensitivity  of  the 
intermediates  such  as  H:  and  CO  to  temperature  variations.  It  was  found  that 
improved  agreements  could  be  obtained  within  the  range  of  the  temperature  uncer¬ 
tainty.  Thus,  modifications  were  again  not  warranted.  In  general,  the  static  reactor 
data  appear  to  be  sensitive  to  the  reactor  temperature  uncertainty. 

Additional  modeling  runs  were  also  conducted  for  the  experimental  oxidation  data 
by  using  the  rate  ofrtsss  and  Tully  (1989)  for  reactions  R157  and  R158.  In  all  cases, 
the  predicted  methanol  profiles  were  found  to  be  steeper  than  the  experimental  data, 
again  indicating  the  need  for  slower  rates  for  these  important  reactions. 
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FIGURE  29  Comparison  between  experimental  (symbols)  static  reactor  oxidation  data  for  ip  =  0.5  and 
$  =  1.0  as  determined  by  Cathonnet  et  al.  (1982)  and  numerical  calculations  (lines)  using  the  present 
scheme  of  Table  I.  (Top:  5.82%  CH,OH.  17.69%  02. 76.49%  N,,  <f>  =  0.5.  Tm  =  823K,p,„  =  0.263  atm, 
isothermal,  isochoric;  Bottom:  5.89%  CHjOH,  8.78%  O,,  85.33%  N2,  <p  =  1.0,  Ta  =  823K, 
/),„  =  0.263  atm.  isothermal,  isochoric:  all  concentrations  on  per  mole  basis). 


8  STUDIES  BASED  ON  SHOCK  TUBE  RESULTS 

Reactor-type  experiments  provide  information  for  the  low-  and  intermediate- 
temperature  regimes,  typically  between  700  K  and  1 200  K.  Shock  tube  studies  comple¬ 
ment  these  experiments  by  operating  at  the  higher  temperature  range  of  1500K  to 
2500  K,  which  is  of  significant  relevance  to  high-temperature  combustion. 

The  present  methanol  mechanism  was  tested  against  the  high  temperature  exper¬ 
imental  results  obtained  from  the  studies  of  Bowman  (1975)  and  Cribb  et  al.  (1984) 
for  oxidation  and  pyrolysis  conditions  respectively.  In  both  investigations  the 
measurements  were  conducted  behind  reflected  shock  waves.  This  configuration  was 
numerically  simulated  as  an  adiabatic,  constant  volume  system,  which  allows  for  both 
temperature  and  pressure  increases  during  the  course  of  the  reaction.  Furthermore, 
the  numerical  results  did  not  require  “time  shifting”  in  order  to  match  the  exper¬ 
imental  data.  The  temperature  uncertainty  for  both  experiments  was  reported  as 
±50K;  the  sensitivity  of  the  predicted  properties  on  this  uncertainty  has  been 
examined  in  all  cases. 

In  Figure  31  the  pyrolysis  experimental  data  of  Cribb  et  al.  (1984)  at  2000K  and 
0.354  atm  are  compared  with  the  numerical  predictions.  In  general,  the  predictions 
were  improved  by  decreasing  the  temperature  to  1950  K,  a  finding  which  is  in 
agreement  with  that  reported  by  Norton  and  Dryer  (1990). 

For  methanol  oxidation  conditions,  the  experimental  data  of  Bowman  (1975)  were 


70 


F.  N.  EGOLFOPOULOS.  D.  X.  DU  AND  C.  K.  LAW 


6 


3 


8 


o 


1 

0 

6 

5 


8 
—  4 


O 

1* 


i 

0 


Time,  sec 


FIGURE  30  Comparison  between  experimental  (symbols)  static  reactor  oxidation  data  for  <t>  =  2.0  and 
<t>  =  4  0  as  determined  by  Cathonnet  el  at.  (1982)  and  numerical  calculations  (lines)  using  the  present 
scheme  of  Table).  (Top.  5.78%  CH,OH.  4.36%  O,.  89.86%  N„  <t>  =  2.0.  Tm  =  823  K.  />,„  =  0.263atm. 
isothermal,  isochoric;  Bottom;  5.78%  CH,OH.  2.13%  O,.  92.09%  N:.  <f>  =  4.0.  Tm  =  823K, 
p,„  =  0.263  atm.  isothermal,  isochoric;  all  concentrations  on  per  mole  basis). 


Time,  psec 

FIGURE  31  Comparison  between  experimental  (symbols)  shock  tube  pyrolysis  data  as  determined  by 
Cribb  el  al  (1984)  and  numerical  calculations  (lines)  using  the  present  scheme  of  Table  I.  (1.0%  CH,OH. 
92  18%  Ne.  1.0%  Ar.  6.0%  H;.  Tm  -  2000K.p,„  =  0.354 atm.  adiabatic,  isochoric;  all  concentrations  on 
per  mole  basis). 
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FIGURE  32  Comparison  between  experimental  shock  tube  oxidation  data  for  <£  =  0.75  as  determined 
by  Bowman  (1975)  and  numerical  calculations  using  the  present  scheme  of  Table  I. 

modeled.  The  parameter  of  interest  is  rma„,  which  is  the  time  interval  between  the 
instants  when  the  reflected  shock  heats  up  the  mixture  and  when  the  product  of  the 
[O]  and  [CO]  concentrations  reaches  a  maximum.  This  parameter  can  be  conveniently 
measured  using  optical  methods  (Bowman,  1975)  and  has  been  determined  under  a 
variety  of  conditions.  All  oxidation  conditions  reported  by  Bowman  were  modeled 
and  close  agreements  between  predicted  and  experimental  values  of  tm,x  were  found 
within  the  stated  ±  50  K  temperature  uncertainty.  The  effect  of  this  uncertainty  was 
greater  in  the  lower  temperature  regime  (1550-1 750  K)  than  for  temperatures  above 
1800  K.  Representative  results  for  mixtures  of  H30H/02/Ar  for  lean  (</>  =  0.75) 
and  rich  (4>  =  6.00)  conditions  are  shown  in  Figures  32  and  33  respectively.  The 
results  indicate  the  discrepancies  at  low  temperatures  as  well  as  the  potential  improve¬ 
ments  obtained  by  varying  the  temperature  within  the  ±  50  K  range.  Similar  behavior 
has  been  observed  for  all  mixtures. 
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FIGURE  33  Comparison  between  experimental  shock  tube  oxidation  data  for  <p  =  6.0  as  determined  by 
Bowman  (1975)  and  numerics'  calculations  using  the  present  scheme  of  Table  I. 
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In  summary,  given  the  temperature  uncertainty  in  the  shock  tube  data,  the  present 
scheme  is  able  to  predict  well  the  experimental  data  in  all  cases  by  simply  varying  the 
temperature  within  the  stated  error  limits.  Because  of  this  agreement,  further  con¬ 
straints  and/or  modifications  of  individual  reactions  are  not  necessary. 


9  CONCLUDING  REMARKS 

In  the  present  investigation,  two  major  objectives  were  accomplished.  The  first  is  the 
experimental  determination  of  stretch-free  laminar  flame  speed  data  for  atmospheric 
methanol/air  mixtures  using  the  counterflow  technique.  These  data  were  obtained  for 
various  mixture  initial  temperatures  and  stoichiometries  ranging  from  very  lean  to 
very  rich,  and  can  be  used  with  confidence  for  the  validation  of  methanol  kinetics  as 
well  as  for  practical  applications.  The  second  objective  is  the  compilation  of  a 
methanol  kinetic  scheme  by  simultaneously  testing  it  against  experimental  data  for 
flame  speeds,  flame  structure,  flow  and  static  reactors,  and  shock  tubes.  The  present 
scheme  is  also  internally  consistent  by  closely  predicting  laminar  flame  speeds  of  mix¬ 
tures  of  CH4  and  the  C, -hydrocarbons  with  air  measured  in  previous  investigations. 

The  modeling  of  various  experimental  configurations  provided  complementary 
information  on  the  kinetics  which  are  essential  for  the  comprehensive  validation  of 
the  mechanism  This  approach  leads  to  the  observation  that  a  rate  slower  than  that 
reported  recently  by  Hess  and  Tully  (1989)  is  needed  for  the  CH,OH  +  OH  -♦ 
Products  reaction  in  order  to  accurately  model  both  flow  and  static  reactor  oxidation 
data.  Such  an  observation  could  not  be  obtained  from  studies  on  laminar  flames 
speeds  alone.  However,  the  possibility  that  certain  reactions  or  combinations  of 
reactions  produce  excess  OH  in  the  fuel  consumption  zone  cannot  be  ruled  out  even 
though  the  present  investigation  could  not  identify  them.  Furthermore,  the  present 
kinetic  scheme  was  also  used  to  identify  sensitivities  and  species  consumption  paths 
in  methanol  flames,  and  to  quantitatively  demonstrate  the  differences  in  for¬ 
maldehyde  and  C:-hydrocarbons  production  between  methanol  and  methane  flames. 

The  problem  of  validating  rates  of  crucial  reactions  with  high  degree  of  confidence 
is.  however,  far  from  being  closed.  In  the  present  investigation,  it  was  shown  that  the 
temperature  uncertainty  for  static  reactors  and  shock  tubes  is  sufficiently  large  that 
the  validity  of  certain  rates  simply  can  not  be  assessed.  Furthermore,  uncertainties 
associated  with  the  rates  of  methanol  reactions  and  branching  ratios,  as  well  as  the 
thermal  decompositions  of  CH;OH  and  CH,0.  leave  substantial  room  for  further 
investigation.  Some  of  these  issues  can  be  clarified  by  studying  systems  at  high 
pressures  at  which  deficiencies  in  our  knowledge  of  thermal  decomposition  reactions 
can  be  rectified. 

Finally,  it  must  also  be  cautioned  that  although  extensive  agreements  were  obtained 
in  the  present  investigation,  the  reported  kinetics  might  not  be  able  to  predict 
methanol  oxidation  and  pyrolysis  properties  when  burning  conditions  deviate  sub¬ 
stantially  from  the  present  range  of  conditions.  The  present  kinetic  scheme,  however, 
does  form  a  useful  starting  point  from  which  further  investigations  can  be  conducted. 
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